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THE DOMAIN OF PARABOLICITY FOR THE MUSKAT PROBLEM 


JOACHIM ESCHER, BOGDAN-VASILE MATIOC, AND CHRISTOPH WALKER 


Abstract. We address the well-posedness of the Muskat problem in a periodic geometry and in 
a setting which allows us to consider general initial and boundary data, gravity effects, as well as 
surface tension effects. In the absence of surface tension we prove that the Rayleigh-Taylor condition 
identifies a domain of parabolicity for the Muskat problem. This property is used to establish the 
well-posedness of the problem. In the presence of surface tension effects the Muskat problem is of 
parabolic type for general initial and boundary data. As a bi-product of our analysis we obtain 
that Dirichlet-Neumann type operators associated with certain diffraction problems are negative 
generators of strongly continuous and analytic semigroups in the scale of small Holder spaces. 
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1. Introduction 

We study the evolution of two vertically superposed (or horizontally adjacent) immiscible layers 
of Newtonian fluids with (possibly) different densities and viscosities in a two-dimensional periodic 
porous medium or Hele-Shaw cell when allowing for both gravity and surface tension effects. The 
pressure on the fixed flat boundary of the lower layer is prescribed and the upper fluid layer is 
assumed to be bounded from above by air at uniform pressure. This leads to a moving boundary 
problem for the interface between the two layers, the interface between the upper layer and the 
air, and the velocity potentials in the two fluid layers. The associated mathematical model is the 
Muskat problem, which was originally proposed in [31] as a model for the encroachment of water 
into an oil sand. It is given in (2.3) below in the absence of surface tension effects and accordingly 
in (7.2) when allowing for surface tension effects. 
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Muskat problems have been studied extensively in the last two decades. Surface tension effects 
are considered in the papers [18, 20, 24, 28] where questions related to well-posedness (for small 
initial data) and stability properties of trivial, that is, circular or flat, and finger-shaped equilibria 
are addressed (see [16] for a classification of equilibria). There is a much larger list of references 
dealing with the Muskat problem without surface tension effects, the methods used in the studies 
being numerous and quite different. The well-posedness property is established in [35,36] by using 
Newton’s iteration method; the references [3,11] use energy estimates (see also [6,12-14] for the 
case of fluids with equal viscosities); [18,20] rely on abstract parabolic theory and the continuous 
maximal regularity due to Da Prato and Grisvard [15]; [33] employs - in the absence of gravity effects 
- methods from complex analysis and a version of the Cauchy-Kowalewski theorem. Existence of 
solutions for nonregular initial data is shown in [5] by means of a fixed point argument. There are 
various interesting phenomena established for fluids with equal viscosities: global existence of strong 
and weak solutions for initial data which are bounded by explicit constants [10,27], existence of initial 
data for which solutions turn over [7-9], or the absence of squirt or splash singularities [13,23,25]. 

An important role in the study of the Muskat problem is played by the Rayleigh-Taylor condition, 
which is a sign restriction on the jump of the gradient of the pressure in normal direction along 
an interface that separates two phases (see (2.7) below for more details) and was originally found 
within the linear theory [32]. In the absence of surface tension effects, the paper [20] was (one of) 
the first in which it was proved that the Muskat problem has, at least for small initial and boundary 
data, a parabolic character provided the Rayleigh-Taylor condition holds. For general initial data 
the well-posedness of the problem in different geometries is also implied by the Rayleigh-Taylor 
condition [3,11,33,35,36]. However, it is worth mentioning that the true character of the problem 
was not revealed in any of the just cited papers. 

In this paper we now prove for arbitrary (sufficiently smooth) initial data that the Muskat problem 
with and without surface tension effects has a parabolic character. More precisely, when neglecting 
surface tension we establish the parabolicity of the problem provided the Rayleigh-Taylor condition 
holds. This enables us to use Da Prato and Grisvard’s abstract parabolic theory, in particular 
continuous maximal regularity, in order to prove the well-posedness of this problem, cf. Theorem 2.1. 
Having two moving interfaces, we actually need to impose the Rayleigh-Taylor at each of them. 
We also show that the Muskat problem with surface tension is parabolic for arbitrary (sufficiently 
smooth) initial and boundary data, the corresponding well-posedness result being stated in Theorem 
7.1. As a bi-product of our analysis we show in Proposition 5.8 and Remark 7.5 that Dirichlet- 
Neumann type operators associated with certain diffraction problems are negative generators of 
strongly continuous and analytic semigroups in the scale of small Holder spaces. 

It is worth to emphasize that the abstract parabolic setting mentioned above appears to be one of 
the few where the Muskat problem can be handled both with or without surface tension effects and 
with general boundary data. In addition, our analysis allows us to handle two fluids with possibly 
different viscosities or densities or even to neglect the effects of gravity (the latter being reasonable, 
for example, when the Hele-Shaw cell is not vertical, but horizontal or also in microfluidic models). 
In many studies these aspects were not taken into account. 

The paper is organized as follows. In Section 2 we present the Muskat problem without surface 
tension and the first main result Theorem 2.1, whose proof requires some preparation. In Section 3 
we first discuss the solvability of a general diffraction problem and recast the Muskat problem as 
a fully nonlinear and nonlocal evolution equation. We then show in Sections 4-6 that the Frechet 
derivative of the operator associated with this evolution problem is an analytic generator. For this 
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we use localization techniques in the spirit of [21], but such that we keep the setting of periodic 
functions. The proof of Theorem 2.1 is then a consequence of this generator result. The well- 
posedness of the Muskat problem with surface tension effects is addressed in Section 7. 


2. The Muskat problem without surface tension effects 

To set the stage we need some notation. In what follows § denotes the unit circle M/(27rZ) 
meaning that functions depending on x € S are 27r-periodic with respect to the real variable x. 
Given m € N and /3 € (0,1), the small Holder space h m+ °(E) stands for the closure of the smooth 
functions C°°(E) in C m+/3 (S). It is well-known that h m+ ^{E) is a true subspace of the classical 
Holder space C' m+/3 (S), cf. e.g. [18], and that C r (§) is densely embedded in h s (E) if r > s > 0. 
Recall that h m (E) = C m (§) for m € N. Similarly, given two functions 4>,ip €: C(S) with x ) 

for all x € § and setting 

H := := {(x,y) : x € S and (j)(x) < y < ip(x)}, (2-1) 

we denote by h m+ ^(Q) the closure of the smooth functions C°°(SY) in C m+/3 (H). As before, C r (H) 
is densely embedded in h s (£l) if r > s provided that £ h s (E). 

Let a € (0,1) and d < 0 be fixed constants and set 

V := £ (h 2+a (§)) 2 : d < f < h}. 

For each pair (f,h) € V we define f2(/) := /) and Q(f,h) according to (2.1). We look for 

a pair of functions (f,h) : [0, To) —» V with To > 0 describing the evolution of the interfaces 
T(/) := [y = f] and F/, := [y = h] that bound two incompressible and immiscible Newtonian 
fluid layers in a porous medium and at constant temperature. At each time instant t € [0, To), the 
domain H(/(t)) is assumed to be occupied by a fluid with density and viscosity /r_, respectively 
H(/(t), h(t)) is the domain occupied by a second fluid with density p + and viscosity p + . Note that 
neither the densities nor the viscosities need to be equal in what follows. We define the velocity 
potentials 

u±:=p± + gp±y, (2.2) 

with g being the Earth’s gravity and p± the fluids’ pressures with densities p±. Our results hold 
also true when neglecting gravity, that is, when g = 0, so we assume g > 0 in the following. The 
velocity fields v± then obey Darcy’s law, see [31], 

k k 

V- =-Vu_ in D(/) and v+ = -Vu+ in Q(f,h), 

p~ p + 

and the incompressibility condition reads 

divu+ = 0 in fi(/) and divu_ = 0 in Q(f,h). 

Here, k > 0 is a constant which stands for the permeability of the porous medium. Assuming the 
pressure on the boundary component T^ := [y = d] to be known, in the absence of surface tension 


4 


J. ESCHER, B.-V. MATIOC, AND CH. WALKER 


effects the Muskat problem is the system of partial differential equations 


Au + 

= 0 

in 

m,h) 

A u- 

= 0 

in 

m), 

d t h 

= — + h' 2 d u u + 

on 

m, 

u + 

= gp+h 

on 

m, 

U- 

= b 

on 

Td, 

U+ — U- 

= g(p+ - p-)f 

on 

rm 

dtf 

= ~kp,± L y/l + f' 2 d v U± 

on 

r (/), 


(2.3a) 


(2.3b) 


governing the evolution of the fluids supplemented with the initial conditions 

/(0) = /o, h(0) = ho. 

We have additionally taken the pressure of the air to be constant zero, we assumed that the inter¬ 
faces between the fluids move along with the fluids, and that the pressure is continuous along the 
interfaces. Given cj) £ C' 1 (§), we have chosen v := (—<//, + (j)' 2 to be the unit outward normal 

vector at the curve [y = <p\. 

The function b = b(t, x ) corresponds to the given pressure at the interface [y = d ] and is assumed 
to belong to the class 

be C([0,T),/i 2+Q (§)), (2.4) 

for some T £ (0, oo]. Our main goal is to study the existence and uniqueness of classical Holder 
solutions to the Muskat problem (2.3), that is, of tuples (f,h,u+,U—) with 

(/, h ) € C([0, To), V) n C 1 ([0, To), (h 1 +a (S)) 2 ) , 

u + (t) € h 2 +a (n(f(t), h(t))), It— (t) £ h 2 +a (tt(f(t))) 

for all t £ [0,Tq) with Tq £ (0, T], and which satisfy the equations of (2.3) pointwise. 


(2.5) 


Given (/o, ho) € V and bo : = 6(0) € h 2+0 
diffraction problem 

A 

A u°_ 

,o 


), Theorem 3.2 and Remark 3.1 below ensure that the 


u + = 
u°_ = 


0 

0 

gp+ho 

bo 


m 

in 

on 

on 


ho), 

mo), 

r (ho), 

r d , 


( 2 . 6 ) 


Um — U°_ 


g(p+ - p-)fo on r(/ 0 ), 

p-d u u° + = p+d u u°_ on T(/ 0 ), 

possesses a unique solution (u^,u((_) £ h 2+a (Q(fo, ho)) X h 2+a (Q(fo)). Letting and p°_ be the 
initial pressures determined, respectively, by u\ and u°_ according to (2.2), we shall show that the 
conditions 

duP- ~ dvP+ <0 on T(/ 0 ), 
d v p° + <0 on T(h 0 ) 

define a regime (for (/o, ho, bo) € V x h 2+a (§>)) where the Muskat problem (2.3) is parabolic. Since 
the air pressure is constant, the condition (2.7) expresses the Rayleigh-Taylor condition imposed 
at each interface as mentioned in the Introduction. To be more precise, we shall prove in Section 
3 that the Muskat problem can be recast as a fully nonlinear abstract evolution equation for the 
interfaces / and h only, that is, 


dt(f,h) = *(t, ( f,h)) 
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which is of parabolic type when (2.7) holds. By parabolicity we mean that the Frechet derivative 
d{f,h)^( 0, (fo,ho)) is the generator of a strongly continuous and analytic semigroup. This property 
is the corner stone in our analysis and, together with the abstract parabolic theory due to Da Prato 
and Grisvard [15,30], it enables us to establish the following well-posedness result for the Muskat 
problem without surface tension effects. 

Theorem 2.1. Let g > 0, (fo,ho) € V, and b be given such that (2.4) holds. Assume that the 
Rayleigh-Taylor conditions (2.7) are satisfied. 

Then, there exist a maximal existence time To := To(fo,ho) € (0,T] and a unique classical 
Holder solution ( f,h,u+,u _) to (2.3) on [0, To). Additionally, the solutions depend continuously on 
the initial data. 

Remarks 2.2. (a) Our analysis discloses that the Muskat problem is backwards parabolic when 
the Rayleigh-Taylor condition holds with reversed inequalities. 

(b) When the fluids have the same viscosities, the last two equations of (2.6) show that the 
Rayleigh-Taylor condition on T(/o), i.e., the first condition in (2.7), is equivalent to p_ > p + (this 
is the case in [ 6 ,12-14])- Hence, in this case the Muskat problem (2.3) is well-posed provided that 

P- > p .|_ and d t ,p° + <0 on T(ho). 

(c) In the particular case d = —1, b( 0) = c € R, and (fo,ho) = (0,1), the conditions (2.7) 
coincide with those found in [18, eq. (2.2)-(2.3)J, that is, 

cu+ , u+ — U- , 

gp+ > - and --- (c- gp + )+g(p+~ p-) < 0. 

P- T+ + P- 

This shows in particular that the set of data (fo,ho,b(0)) for which the Rayleigh-Taylor condi¬ 
tions (2.7) are satisfied is not empty. 

(d) If gravity is neglected, that is, if g = 0, the Rayleigh-Taylor conditions (2.7) are equivalent 
to 

(p+ ~ fi-)d v P+ < 0 on r(/ 0 ), 
d u p+ <0 on T(ho). 

As for Remarks 2.2 (d) we point out that if bo is zero or a negative function, then (2.8) cannot 
be satisfied. Indeed, if bo is the zero function, then both p° + and p°_ are identically zero. If bo 
is a negative function, then p°_ is also negative since otherwise there exists xq € § such that 
p°_(xo, f(xo)) = max^_ p°_ > 0, and so d u p°_ > 0 at (xo,f(xo)) by Hopf’s lemma. But p° + is 
harmonic as well and not constant, therefore p+(xo, f(xo)) = max^ p { )_ > 0, so that < 0 at 

(xo, f{x o)) in contradiction to the last equation of (2.6). Hence, p°_ is negative implying p® + < 0 on 
r(/o). By Hopf’s lemma we find d v p® + > 0 on T(/io), and (2.8) is again not satisfied. 

Lastly, if bo is a positive function, the previous arguments show that (2.8) is equivalent to 

(T+ ~ P-)dvP+ <0 on r(/ 0 ). 

As p° + attains its positive maximum on r(/o), we have that d u p q” is negative at least at one point 
on this interface implying that 

p+ < n— (2.9) 

It is easy to see that if bo is positive and constant and if also fo and ho are constant functions, then 
(2.9) is equivalent to (2.7). Thus, if bo is positive and constant we find, as in [33], that the Muskat 
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problem is well-posed for small initial data - that is, initial data close to constants in (/i 2+ “(§)) 2 - 
when the more viscous fluid expands into the less viscous one. 

The proof of Theorem 2.1 is postponed to the end of Section 6 as it requires several preparatory 
results that will be given in the subsequent sections. 


3. The evolution equation 

In order to solve problem (2.3) we re-write it as an abstract evolution equation on the unit circle. 
To do so we first transform system (2.3a) into a system of equations on fixed domains by using the 
unknown functions ( f,h ). Let := § X (—1,0), 0+ := § x (0,1), and define for each (/, h) € V 
the mappings </>/ : -A 0(/) and 4>(f,h) • D(/, h) by setting 

4>f(x,y) ■= ( x,-dy + (1 + y)f(x)) and (f>^ h) (x,y) := (x,yh(x) + (1 - y)f{x)), 

respectively. One easily checks that cj)f and are diffeomorphisms for all (f,h) € V. Each pair 

(/, h) € V induces linear uniformly elliptic operators 

A(f) : /i 2+ “(fl_) -A- h"(0_), V- ha A(u_ o fj 1 ) o (f>f, 

A{f, h) : h 2+a (n + ) -A h a (n + ), v + ha A(u+ o ^ h) ) ° 

which depend, as bounded operators, real-analytically on / and h (see the formulae in the Appen¬ 
dix). Denote by tro the trace operator with respect to To := § x {0}. We associate with problem 
(2.3a) trace operators on To, 

B(f)v- ■= kpZ 1 tr 0 ((V(u_ o 1 )|( 1)) o (/>/), v _ € /i 2+ “(D_), 

B(f,h)v+ ■■= /c/r+ 1 tro((V(u + o 0-^)1 (-/', 1)) o <£ (/j/l) ), v + G h 2+a (Q + ), 

which, seen as bounded operators into h l+a { §), depend real-analytically on / and h as well. Lastly, 
we define a boundary operator on Ti, where T±i := § x {±1}. Given (/, h) G V, we set 

&i(f, h ) v + ■= k V+ l tri«V(u+ o cj)-^ h) )\(-ti, 1)) o 0 (M} ), v + € h 2+a (n + ), 

where trj-i is the trace operator with respect to T±i. 

Remark 3.1. Given ( f,h ) € V, the mappings 

[u- i —y U— o (j)f\ : h 2+a (£!(/)) -A /i 2+ “(D_), 

[u + Ati + o (j) {fA) ] : h 2+a (n(f, h)) -A h 2 + "(D+), 


are isomorphisms. 

Proof. See, for instance, the proof of [19, Lemma 1.2]. 


□ 
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In view of Remark 3.1, it follows that (/, h, u + , uJ) is a solution to (2.3) if and only if (/, h, v + ,v_) 
with v+ := u+ o and v - := u - ° is a classical Holder solution to 


(3.1a) 


•A{f, h)v + 

= 0 

in n + , 


= 0 

in ff_, 

~B(f) v- 

= 0 

on T 0 , 

V + — V- 

= g{p+-p~)f 

on T 0 , 

V+ 

= gp+h 

on Ti, 

V- 

= b 

on T_i 


with 


and 


d t h = -Bi(f,h)v + on Ti, 
dtf = -B{f)v- on r 0 , 


/( 0 ) = /o, h( 0 ) = ho. 

The notion of classical Holder solution to (3.1) is defined analogously to that for problem 


(3.1b) 

(3.1c) 

(2.3). 


A diffraction problem in Holder spaces. The system (3.1a) is an elliptic diffraction (or 
transmission) problem, problems of this type being highly relevant in many physical situations such 
as the study of multiphase dynamics. However, citable references on this topic are sparse. The main 
goal in this part is to establish the following result on the existence, uniqueness, and real-analytic 
dependence of solutions to (3.1a.) on given (f,h) G V and b G h 2+a (§). 

Theorem 3.2. Given ( f,h ) G V and b G h 2+a (E>), there exists a unique solution 

(«+.«-) := («+(/, M),«-(/,M)) G h 2+Q (n + ) X h 2+Q (n.) 

to the diffraction problem (3.1a). Moreover, it holds that 

[(/, h, b) ha (v+(f,h,b),v-(f,h,b))] € C“(V x h 2+a (S),h 2+a (Q+) x h 2+a (n_)). 

In the Holder setting considered herein, problem (3.1a) can be accessed by using the celebrated 
Agmon-Douglis-Nirenberg estimates on solutions to elliptic boundary value systems presented in [1]. 
To prove Theorem 3.2 we first consider a particular boundary value problem for a linear elliptic 
system with coupled boundary conditions for which we establish the existence and uniqueness of 
solutions in the natural Holder spaces. This is the context of the next proposition. 

In the following d\ := d x , 82 := d y and we identify both boundaries of f2+ with the unit circle S. 


Proposition 3.3. Let := dij+b^di+c^ be uniformly elliptic operators with coefficients a[j\ 
b[ k \ c( fc ) € C Q (H + ) for i,j, k = 1, 2 and such that c^ < 0. Additionally, let Bj. := (3^ tro <9* + 7 ^ 
be two boundary operators such that fd[ k \ 7 ^ G C 1 +Q (S), f}^ > 0, and 7 ^ < 0, k = 1,2. Given 
Fi,F 2 G C“(f2+), tp 1 G C 1+ “(§), and p 2 -,P'Si l Pi G C' 2 +a (S), the boundary value problem 


and 



( Ciwi = Fi 

in fi + , 



\ C 2 W 2 = F 2 

in H + , 


B\w\ + B 2 w 2 

= pi on T 0 , 

f w\ = 

= p 3 onTi 

W\ — W 2 

= ip 2 on T 0 , 

1 W2 = 

= onT 1 


possesses a unique solution {wi,wf) G (C' 2 + “(H_ ) _))" 


(3.2a) 


(3.2b) 
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Proof. We discuss uniqueness first. To this end, let (w\,w 2 ) be a solution to (3.2) with right-hand 
sides replaced all by zero. If max^ w\ > 0 (the case min^ + w± < 0 is similar), the weak elliptic 
maximum principle ensures that 

inaxioi = max w 2 = rt>i(x,0) = w 2 (x, 0 ) 

Q_l_ Q_|_ 


for some x G §. Applying Hopf’s lemma at (x,0) for both w\ and w 2 yields d y wi(x,0 ) < 0 and 
d y w 2 ( 2 , 0 ) < 0. This is in contradiction to the equation B\W\ + B 2 W 2 = 0 on To- Hence, (w\,w 2 ) 
has to be the zero solution, and therefore (3.2) has at most one solution (wi,ii] 2 ) G (C' 2+Q (fl+))T 
For the existence part, we consider a family of operators {71 '■= (T) 1 ,..., 7i( 6 )} rg [o ) i] C £(X, Y), 
with X := (C 2+ "(D + )) 2 , ¥ := (C"(D+)) 2 x C 1+ “(S) x (C 2+a (§)) 3 , and 


T t (wi,W 2 ) 


( (1 — t)L\W\ + tAw\ ^ 

(1 - t)C 2 w 2 + tAw 2 

(1 - t)[BiW\ + B2W2) + Ttr 0 d y (wi + w 2 ) 
tr 0 (wi - w 2 ) 
tri w 1 

V tri w 2 / 


for r G [0,1] and ( 101 , 102 ) G (C' 2+ “(f2+))T We observe that [r eA 71] G C([0,1], £(X, ¥)). Moreover, 
by considering two suitable Dirichlet problems the invertibility of 71 can be reduced to the solvability 
of the equation 71 ( 21 , 22 ) = (0,0, p, 0, 0,0,0) for arbitrary p G C 1+a (§). Indeed, letting ( 21 , 22 ) 
denote the solution to the equation Ti(zi,z 2 ) = (0, 0, ip, 0,0, 0,0) for ip := p\ — tro d y {w\ + 102 ), and 
setting 

:= (A, tr 0 , tr 1 ) _ 1 (l ? i, (/? 2 , ^ 3 ) and w 2 ■= (A, tr 0 , tri) _1 (T 2 ,0, tp 4 ), 

it is easy to see that (wi,w 2 ) := (21 + w\, z 2 + w 2 ) G X satisfies Ti{w\,w 2 ) = {F\,F 2 , pi,p 2 ,P 3 , <£> 4 ). 
The solution to T\{z\, z 2 ) = (0,0, p, 0, 0,0, 0) is 21 = z 2 '■= (A, tr 0 d y , tri) _ 1 (0, p/2, 0) € C 2 +a (PL + ). 
Hence, we have shown that 71 is invertible. If we find a constant C > 0, such that 


||(wh,io 2 )||x < C\\T t (wi, io 2 )||y for all r G [0,1] and (ioi,io 2 ) G X, (3.3) 

then by the method of continuity, cf. e.g. [26], we conclude that 7 o is an isomorphism, which is the 
claim of the proposition. 

We are left to establish (3.3). To this end, we show that 71 corresponds to a uniformly elliptic 
system that satisfies the Complementing Condition in the sense of [1] on both boundary components 
To and Ti. Because the equations in are decoupled, it is easy to see that defines 

a uniformly elliptic system for each r G [0,1]. Additionally, the boundary conditions defined by 
(71,7)?) are of Dirichlet type, and therefore the Complementing Condition on Ti is straightforward. 
To verify the Complementing Condition on To we modify the operators 71*, 1 < i < 4 as follows: 
we identify the principle parts Tt' 1 , 1 < i < 4, freeze their coefficients at an arbitrary P G To, and 
replace (di,d 2 ) by (£, —id/) with 0 / ( G M. Doing this, we arrive at the initial value problem 

v'[ — iA^v^ — A^v 1 = 0, for t > 0, 

v'{ — iA^v ' 2 — A^v 2 = 0 , for t > 0 , 

< 'Ui(O) = u 2 (0), 

*[((! -t)$\P) + r K(°) + {i 1 - t )^ 2 2) ( p ) + T ) V 2( Q )} 

-(1 -r)C(/3j 1) (P)u 1 (0) + /3f ) (P)u 2 (0)) = 0, 


(3.4) 
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where 


:= - 


2(1 -r)ag(P)C 


Ak) ((1 ~r)a^(P) + T)e 


,( fc )/ 


k = 1 , 2 . 


(1 — r)a^(P) + t " (1 - r)a^(-P) + t 

The Complementing Condition is satisfied if and only if the only bounded solution (v\, ry) of (3.4) 
is the zero solution. It is readily seen that 

v k (t) = k = 1 , 2 , 


with := /2 and 6 ^ := \J — (H*^) 2 /4 > 0. The boundedness of V\,V 2 entails that 

7 ^ = 7 ^ = 0. Moreover, the equation ui(0) = 02 ( 0 ) implies that yj 1 ^ = 7 ®. Finally, assuming 
yj 1 ) 7 ^ 0, we hnd from the last equation of (3.4) that necessarily 

4 1} ((! - r)4 1] (P) + r) + 4 2) ((1 - t)/4 2) (P) + r) = 0. 

However, this last equation cannot hold true as <51, ' and /3 2 (P),k = 1,2, are positive constants. 
Hence, the Complementing Condition is also satisfied on To- 

We may use now Theorem 9.3 and argue similarly as in the subsequent Remark 2 in [1] to 
conclude, together with the uniqueness result established at the beginning of the proof, that there 
exists a constant C > 0 such that the estimate (3.3) holds. This completes the proof. □ 


Using Proposition 3.3, we obtain the unique solvability of certain diffraction problems within the 
natural Holder spaces. 

Corollary 3.4. Let C± := flydy + bfdi + cF* 1 be uniformly elliptic operators with coefficients 
dy , bf , G C a (Ll±) and such that < 0 for i,j = 1,2. Moreover, let B± := flf 1 try d t + y 1 * 1 be 
boundary operators such that G C 1+Q (§), (3 2 > 0, and y 1 * 1 < 0. Then, given F± G C a (Q±), 

ip 1 G C' 1+ “(§), and 992 , £3, £4 £ C 2+ "(§), the diffraction problem 


(3.5) 


possesses a unique solution (u+,u_) G C' 2 “ l ~ a (fl+) x C 2+ “(fl_). In particular, there exists a constant 
C > 0 such that 


C + v+ 

= F+ 

in H + , 

£_u_ 

= F_ 

in H_, 

-B^v- 

= £1 

on T 0 , 

V + — V- 

= £2 

on T 0 , 

V+ 

= £3 

on T1, 

V- 

= £4 

on T_i 


llwlb+o + ||t>—||2+a < C f ||F+|| a + \\F_\\ a + 


A\\i+ a + ^ I 

i =2 


f>i\\2+a 


(3.6) 


Proof. The mapping <f> : Q_)_ —>• dehned by = —y, y G f2+ is a smooth diffeomorphism. 

Therefore, if (v+,V-) solves the system (3.5), then w\ := v + and W 2 = V- o </>, is a solution to (3.2) 
with £1 := £ + , £2 := [w (£_(ro o </> -1 )) o </>], B\ := £>+, B 2 '■= [re >-)• —B-(w o ^ _1 )], F\ := F + , 
and F 2 := £_ o f. The desired claim follows now directly from Proposition 3.3. □ 


We are now in the position to prove Theorem 3.2. 
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Proof of Theorem 3.2. Given (f,h) G V, it follows from Corollary 3.4 that the mapping 


(v+,v-) 


/ A(f,h)v+ \ 
•A(f)v- 

B{f,h)v+-B(f)v- 
tr 0 (u+ - v-) 
tri v + , 

\ tr_i V- / 


(3.7) 


defines an isomorphism between C 2+a {PL + ) x C 2+ "(fl_) and x C a (Q-) x C' 1+a (S) X 

(C 2+a (S)) 3 . Since (f,h) G V, a density argument shows that the operator defined by (3.7) is a 
isomorphism also when acting between the corresponding small Holder spaces h 2+ "(H + ) x /i 2+ “(H_) 
and h a (£l+) x /i a (fl_) X /i 1+ "(§) x (h 2+a ( S)) 3 . Because the differential operators and the right-hand 
sides of the equations of (3.1a) depend in a real-analytic way on (/, h, b) G V x h 2+a (S), the claim 
of Theorem 3.2 is now obvious. □ 


The evolution equation. With the identification T* = § for i G { — 1,0,1} and by using The¬ 
orem 3.2, the problem (3.1a)-(3.1b) can now be reformulated as an abstract fully nonlinear and 
nonlocal evolution equation 

= (3.8) 

where <f> : [0, T) X V —» ( h 1+a (S )) 2 is the operator 4> = (<f>i,<f> 2 ) defined by 

* 1 (t,(f,h))--=-B(f)v-(f,h,Kt)), 

:=-Br (f,h)v+(f,h, b(t)), ( j 


with (v+,v-) denoting the solution operator introduced in Theorem 3.2. We note that 

$eC([0,T)xV,(/i 1+a (§)) 2 ) and $(t, ■) G C“(V, (/i 1+a (S)) 2 ) fortG[0,T). (3.10) 


Let (/o, ho) G V and set 

b 0 := 6(0) G h 2+a (S). 

Our aim is to apply the existence result [30, Theorem 8.4.1] to (3.8) for which we need to show that 

- d {W m(fo,h 0 )) €n{(h 2+a (S)) 2 ,(h 1+a (S)) 2 ), (3.11) 


that is, <9 (j m 4>(0, (/o,^o)) seen as an unbounded operator in (/i 1+ "(§)) 2 with domain of definition 
(h 2+Q (§)) 2 is the generator of a strongly continuous analytic semigroup in £((/i 1+ "(§)) 2 ). This 
generator property will be established in the sections to follow for (fo,ho) G V and bo G h 2+a (E>) 
for which (2.7) is satisfied. 

Given (/*,h*) G V, the operator 5(/,ft)^(0, (/*, h*)) G C[(h 2+a (S)) 2 , (h 1+a (§)) 2 ) can be written 
in matrix form 


(/•>*.)) 


/<9/4>i(0, (/*, h*)) d h $i (0, (/*, K))\ 

\ v «9 / T» 2 (0, (/*, h*)) d h $ 2 (0, (/*, h *))) 1 


where, according to the definition (3.9), we have 

df* i(0, (f,,h,))[f\ = -d f BUMv-U*,h,,b 0 )-BU*)dfV-(f,,h,,b 0 )[f\, 

dh.® i(0, (f*,h*))[h] = -B(f*)d h v-(f*,h*,bo)[h], 

dh$ 2 ( 0 , (/*, /i*))[/i] = -d h Bi(f*,h*)[h]v+(f*,h*,bo) - #i(/*, h*)d h v+(f*, h*, b 0 )[h], 


(3.12) 
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for (/, h) G (/i 2+ "(§)) 2 . Additionally, (w+[f],W-[f]) := (dfV + (f*,h*,bo)[f],dfV-(f*,h*,b Q )[f]) is 
the solution to the diffraction problem 


A(f*,h*)w+{f] 

= -d f A{f*,K)[f]v% 

in 

A(f*)w-[f] 

= -d s A{f*)[f]v*_ 

in H_, 


= -d f B(U,K)\J]v% + d f B{U)\J}v*_ 

on T 0 , 

™+lf] ~ w-[f] 

= g(p+-p~)f 

on T 0 , 

w+[f\ 

= 0 

on Ti, 

«>-[/] 

= 0 

on T_i 


and (W+[/i], W_[/t]) := (c^ff^/*, /i*, 6 o)[^]j <9fcT* (/*, h*, bo)[h\) solves the diffraction problem 

A(f*,h*)W+[h] = -dhA(f*,h*)[h]v+ 

A(f*)W-[h] = 0 

B{U,K)W+[h]-B{U)W-[h] = -d h Btf*,K)[h}v% 

' W+[h]-W-[h] = 0 

W+[h} = gp+h 
W-[h] = 0 

In (3.13) and (3.14) we set 

(*>+,«*) = {v+,v-)(f*,h*,b 0 ). 

According to [2, Theorem 1.1.6.1 and Remark 1.1.6.2], (3.11) is satisfied provided the diagonal 
operators satisfy 

- 9/^(0, (/ 0 ,V>) € W(/» 2+a (S),/i 1+a (S)), (3.16) 

- d h Mo,(f 0 ,h 0 )) en(h 2 +a (§),h 1+a (§)), (3.17) 

and provided the following property holds for the off-diagonal operator: for each e > 0 there exists 
A'o = K 0 (e) > 0 such that 

||9^!(0, (/o, ho))[h}\\i +a < e\\h\\ 2+a + K 0 \\h\\ 1+a for all h G h 2+a ( S). (3.18) 

So, to establish (3.11) various computations are needed. In Section 4 we first prove (3.18) based on 
Schauder estimates for diffraction problems as presented in Corollary 3.4. The proofs of (3.16) and 
(3.17) respectively, are given in Sections 5 (see Theorem 5.7) and 6 (see Theorem 6.5) for which 
we use localization techniques in the spirit of [21] (see also [17,22,37]). However, our localization 
techniques are quite different from those therein as we do not consider problems in the halfplane, 
and our results are sharper, e.g. see Theorem 6.5 and [37, Theorem 14]. 


in 

n + , 

in 

n_, 

on 

r 0 , 

on 

r 0 , 

on 

n, 

on 

r_i 


(3.15) 


4. An off-diagonal operator 

The main goal of this section is to establish the property (3.18). This is a consequence of the 
following lemma where (3.18) is established for general (/*,/i*) G V. 


Lemma 4.1. Let (/*, /i*) G V. Given s G (0,1), there exists Kq = Kq(e) > 0 such that 
H^rtO, (/*, K))[h\\\ 1+a < e\\h\\ 2+a + K 0 \\h\\ 1+a for all h G h 2 +a ( S). 


(4.1) 
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Proof. The proof is based on estimates for boundary value problems for elliptic systems, cf. (3.6) 
and [1], Let (/„, h n , b n ) be a sequence in (C°°(§))’ which converges towards (/*, fi*, b o) in (/i 2+a (§)) 
and such that (/„, h n ) G V for all n. We then have 

||d h $i(0, (f*,h*))[h] ||i +Q , < C\\ tr 0 (Vd h V-(f*, K,b 0 )[h\)\\ 1+a 


< C || tr 0 V (d h v-(f n , h n , b n )[h ] - d h V-(f*,K,b 0 )[h]) ||i +Q + C\\ tr o(Wd h v-(f n , K, b n )[h])\\i +a 
— C\\dhV-(f n , h n , b n ) — dhV-(f*, h *, ^o)||£^ 2 +a(g) ll^lb+a 

+ C\\ti 0 V(d h V-(f n ,h n ,b n )[h])\\ 1+a =: h + I 2 (4.2) 

for all h G h 2+a ( S). Let e > 0 be given. In view of Theorem 3.2, we can choose n large enough to 


guarantee that 


h<- 
“ 2 


2+a- 


(4.3) 


We are now left to estimate the term I 2 for a fixed n such that (4.3) holds. To this end, we 
let (W^[h\,Wf[h\) := (d h v + (f n ,h n ,b n )[h\,d h V-(f n ,h n ,b n )[h\) denote the solution of (3.14) when 
replacing (/*,fi*,6 0 ) by (f n ,h n ,b n ) and v* + by vf := v+(f n , h n , b n ). Because (f n ,h n ,b n ) is smooth, 
it follows from [1, Theorem 9.3], by arguing as in the proof of Theorem 3.2, that vf is a smooth 
function up to the boundary of f }+. In particular, vf G /i 3+ "(fl + ). We now split the solution 
(W%[h],W”[h]) as 

(W$[h ], W?[h ]) = (W+, Wl) + (W+, Wl), (4.4) 

where (W\, W}_) is the solution to 


(4.5) 


respectively solves 


A(f n ,h n )Wl 

= —d h A(f n , h n )[h\vf 

in 

A (f n )Wl 

= 0 

in fl_, 

r i ~ B(f n )Wl 

= 0 

on T 0 , 

w\-wi 

= 0 

on T 0 , 

wi 

= 0 

on Ti, 

wi 

= 0 

on T_i 

A(f n ,h n )Wf 

A(f n )Wf 

= 0 

in 

= 0 

in II_, 

T l ~ 

= -d h B{f n ,h n )[h}vl 

on r 0 , 

Wl - Wf 

= 0 

on r 0 , 

W 2 

= 9P+h 

on Tr, 

wf 

= 0 

on r_i 


(4.6) 


In order to estimate Wf we note that 

d h B{f n ,h n )[h} = -kfrf 1 j- + ^ n ^ tr 0 d y , 

\^n Jn) 

and therefore, the right-hand side of the third equation of (4.6) belongs to h 2+Q (§). In virtue 
of [1, Theorem 9.3], there exists a constant C > 0 (which depends on the previously fixed n) such 
that on the subdomain (l/2)fi_ := § x (—1/2,0) of we have 

||W2||^ )n - < C(\\d h B(f n ,h n )[h]vl\\ 1+a + \\Wl\\o + \\W*\\ 0 ) 


(4.7) 
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for all h G h 2 +a (S). Let a' € (0, ct) be fixed. Recalling (3.6), we have 

Iloilo + II^Ho < IIW^II 2 +a' + \\W*\\ 2+a' < C(\\d h B(f n ,h n )[h]vl\\ 1+a , + \\h\\ 2+a ,), 
and together with (4.7) we end up with 

\\W 2 _\\^~ < C(\\d h B(f n ,h n )[h]vl\\ 1+a + \\h\\ 2+a ,) < C\\h\\ 2+a/ . (4.8) 

Using the following interpolation property of the small Holder spaces (e.g. see [30]) 

(h r {S), h s (§)) e = fcMH-^s) fo r 0 G (0,1) and (1 - 6 )r + 6 s£ N, (4.9) 

where (•, -)g = (•,•)() is the real interpolation functor introduced by Da Prato and Grisvard [15], 
we infer from (4.8) and Young’s inequality that there exists a constant K(e ) such that 

||tr 0 VW2|| 1+a < \\W 2 _\\^l )n - < ^\\h\\ 2+a + K(e)\\h\\ l+a (4.10) 

for all h G h 2+a (§). 

Thanks to (4.10), we are left to prove a similar estimate for || tro VWl||i+ a . To this end, we 
choose 5 G (0,1) and a function x '■= XS € C' oo ([0,1[) such that 0 < x < 1, y = 1 on (0,5/8), and 
X = 0 on [5/4,1], Note that (3.6) and the regularity A G C“(V,£(/i 2+ “(S), /i 1+Q (§))) implies that, 
for a fixed a' G (0, a), 

llWlHa+a' + ||Wj|| 2+a , < C\\h\\ 2+a , (4.11) 

for all h G h 2 +a (S). Using the more compact notation A(f n , h n ) = d, t j + b 2 d 2 , the pair (W+x> LUl) 
solves the system (4.5), but with the first equation replaced by the elliptic equation 

A(fn, h n ){W l + x) = - xd h A(f n , h n )[h]v n + + 2 audAVlx! + a 22 {d 2 W\^ + W^x") + Wjx' 
in n_|_. We now use (3.6) and (4.11) to get 

WWih+a <C^\xd h A{f n ,h n )[h\vl\\ a + \\{2a l 2 d l Wl + b 2 Wl + a 22 d 2 Wl)x'+ a 22 Wlx''\\o) 
<C\\xd h A(f n ,h n )[h\vl\\ a + C( 6 )\\Wl\\ 1+a 
<C\\yh" X \\ a + C (5)\\h\\ 2+a t 

<C\\yx\\o\\hh + a + C( 6 )\\h\\ 2+a '. (4.12) 

Choosing 5 such that C||yxl|o < s/4, we infer from (4.12) that there exists a constant K(s) such 
that 

||tr 0 VWi:||i +a <||Wi:|| 2 + a < (e/4)||/i|| 2+a +iir(e)||/i|| 1+a (4.13) 

for all h G h 2+a (§). Gathering now (4.2), (4.3), (4.4), (4.10), and (4.13) we have established the 
desired estimate (4.1). □ 


5. The first diagonal operator 

In this section we prove that <9f<Fi(0, (/q, ho)) is the generator of a strongly continuous and 
analytic semigroup, hence (3.16), when (fo,ho) and ^0 are such that the first inequality in (2.7) 
holds. This is stated in Theorem 5.7. We start in Lemma 5.1 by identifying the “leading order part” 
of 9/-4>i(0, (/*,/i*)) for a general (/*,/i*) G V. This reduces our task (see (5.4) and [2, 
Theorem 1.1.3.1 (ii)]) to showing the generator property merely for df^\{fo, ho), and thus allows us 
to neglect several “lower order terms” in the quite involved computations to follow. Following this 
step, we locally approximate the principal part dy (0, (/*,/i*)) by Fourier multipliers for (/*,/i*) 
sufficiently close to (/q, ho) and possessing additional regularity, cf. Theorem 5.5. These multipliers 
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are then shown to be generators of strongly continuous analytic semigroups, where the constants 
in the resolvent estimates are uniform with respect to certain variables, cf. Lemma 5.6. This 
uniformity property is essential when establishing the desired Theorem 5.7 by means of Lemma 5.2 
and a continuity argument. 


Lemma 5.1. Given € V, let df$>f(f*,h*) € £(/i 2+ “(§), /i 1+ “(§)) denote the operator 

defined by 


d f n(.f*,h*)[f] '■= — 


k /2fitr 0 d v vl 


-tr 0 d x v*_)f'-B(f,)w*_[f], /a 2+ “(S), (5.1) 


li- V f*~d 

with (v+,vf) defined in (3.15) and where (wT [/], uU [/]) denotes the solution to the problem 

AUf*,h*)wl[f] = ^° dyV r + f" 




h * - f, 
tr 0 d v v* 


f* — d 


V ~f" 


w+[f] ~ w-lf) = s{p+ ~ P—)f 
wlif] = 0 
W-[f] = 0 


with 


Aq (/*) d xx —t—- -,d X y + 77-*4.q(/*j h*) d xx — — - ~rd X y + 


in 

n + , 

in 


on 

r 0 , 

on 

r 0 , 

on 

Ti, 

on 

r _ l5 

f? + 1 


(5.2) 


f*-d xy (/*-dy 

d f B*(f*)[f] := —(T^tro^-tro d x )f, d f B"(f*,h*)[f] ■= — ( 

U— \ /* — CL / U -u \ 


:d v 


' {K-UY yy ' 

k ( 2 fi 


P-\- ^ h* f* 


tr 0 d y — tr 0 d x ^j /'. 

(5.3) 


d 

Then, given e € (0,1), there exists K\ = Lfi(e) > 0 such that 

II(0, (/* A)) [/] - d f $l (/*, /»*)[/] II 1+a < £||/|| 2+ a T^ll/lk+a for all f £ h 2 +“(S). (5.4) 


Moreover, d f <$>\ € (V, £(/i 2+q (§),/ i 1+ “(§))). 

Proof. The regularity property follows by using similar arguments as in the proof of Theorem 3.2. 
In order to prove (5.4), let £ € (0,1) be given and a! € (0, a) be fixed. We choose 5 £ (0,1) and 
a cut-off function X := XS € C°°([— 1,1]) such that 0 < x < 1, x = 1 f° r |y| < 5/2, and x = 0 for 
|y| > 5. Using the same notation as in Section 3 (see (3.13)) we have from (3.12), (5.1), and the 
formula for djB(f *) in the Appendix A 

l|3/^(0, (/.,/*.))[/] -df*Uf*,h*)[f]\\i+c l <C(||/|| 1+a + ||S(/*)(n;_[/]-^[/])|| 1+Q ) 

<C(\\f\\l+a + || tr 0 V (w— [/] - «£[/]) ||i +a ) 
<C(||/||i + „ + Hx(»-[/] - ! 0 l[/])ll? + -J, (5.5) 

where C is independent of 5. We now observe that 


(«+,«-) : = i.wl[f],wl[f]) - (w + [f],w-[f}) 
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solves a diffraction problem of the form 

A%(f*,h*)u + = a^f + af f' + ya\ ((1 — y)d y v* + — tr 0 d y v* + )f" 

+y(bod xy w + [f] + bfdyyW + [f]) + b^d y w + [f] in Q + , 

Aq {f*)u- = a^f + a^f + yai + ((1 + y)d y v*_ - ti 0 dyV*_)f" 


+y(b 0 d xy w-[f} + b 1 dy y w-[f ]) + b 2 dyW-[f] in S4_, 

B(f*,h*)u+- = cf on r 0 , 

u + — U- = 0 on To, 

u + = 0 on Ti, 

U- = 0 on T_i, 


with functions af,bf € h a (£l±), 0 < i < 2, and c € h 1+a (§). Therefore, (u*,u*) := x( u +i u -) is 
the solution to 


Al(f*,h*)u\ 

AqU*) u - 

< B{U,K)u\-B{U)u x _ 

u x — uf 


xA%(f*, K)u + - -j^jrx'dxU+ + (2x?d y u+ + x"u+) 
xAUf*)u- ~ jy^x'dxu- + (‘ 2 x'dyU- + x"u~) 

cf 


0 

0 

0 


in Q + , 
in fi_, 
on T 0 , 
on T 0 , 
on Ti, 
on T_i 


We can estimate the solutions to the systems (5.2) and (5.6) by using (3.6) and obtain that 


(5.6) 


u x -fc <c(\\xAUf*,h*)u+\\%+ + \\xA5(f*)u-\\%- + \\f\\i +a )+C(8)(\\u + \\ 2+al + ||u_|| 2+a ,) 
<c(||x((l - y)d y v* + - tr 0 d y v* + )f"\\a + + Ilx((l + y) d v v - - tr o d y v*_)f "||?" 

+ \\xy{b^d xy w + [f}+b+d yy w + [f])\\2 + + \\xy{bo d xy w_[f] + b^d yy w_[f})\\^ 

+ C( 5 )\\fh+a> 


< c '(llx(( 1 -y)dyV* + -tr 0 d y vX)\\^ + + Wx{(l + y)dyV*_-ti 0 dyV*_)\\Q 

+ Ilx 2 /ll? + + llxyll?') ll/112+a + c(S)\\f\\ 2+a , 

Recalling the definition of y = xs and choosing 5 > 0 such that 

llx((l - y)dyV* + -tr 0 d y v* + )\\Q + + \\x((l + y)d y v*_ - tr 0 d y v*_)\\^- + \\xy\\o + + 11x2/11?' < 
the desired estimate (5.4) follows from the interpolation property (4.9). □ 


Let (fo,ho) and bo be such that the first inequality of (2.7) holds. We are now left to show 
that the principal part df&J(fo,ho) of c//<I>i(0, (/o, ho)) is the generator of a strongly continuous 
and analytic semigroup in £(h 1+ "(§)). In view of the classical result [2, Theorem 1.1.2.2], one has 
—c//4>i (/o, ho) € L/(/i 2+ “(S), h 1+Q (§)) if and only if there exist constants K\ > 1 and > 0 such 
that —df&f(fo,h 0 ) € %(/i 2+a (§), h 1+a (S), Ki,wi), that is, 

wi - h 0 ) € Isom(h 2 +“(§), h 1+a (8)) 
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and 


< « for -11 ReA > M and 0 / / £ fc 2 +“(S). 


|A| • ll/lll+a + ll/lb+a 

To this end, for cr > 0, let S a denote the set consisting of those (/*, h*) G V satisfying the inequalities 

(1) a < min{/* - d, h* - /*}, \\f*h + IIM 2 < cr -1 , 

, \ tr o d y v*_ tr 0 d y v* + 

(2) 9 0_^ + )>n +7 _ 7 -- 7 _±, 

(3) || tr 0 dyV+\\o + || tr 0 «9 y i;l||o < o'” 1 , 

where {v’f^v'f) is defined in (3.15). Since the functions (/o, ho) and b q are chosen such that the hrst 
inequality in (2.7) is satished, we may choose a such that (fo,ho) G S a . 


Lemma 5.2. Let a > 0 be such that (/o, ho) G S a . Assume that there exists a constant := ki(cj) 
and for each (/*,h*) G (h 3+ “(§)) n<S , 0 -nB(/ l 2 +a(g ))2 ((/o, ho), cr) there exists a further constant uj\ > 0 
with the property that 

-d f ^(f„K) en(h 2 +a (S),h 1 +a ( S),«i,wi). 

Then 

-d f *Ufo, h 0 ) G n{h 2+a (§), h 1+a (8)). 

Proof. Using the regularity assertions in Theorem 3.2 and Lemma 5.1 together with the density of 
/i 3+ “(§) in h 2+a (S), we find (/*,h*)G (/i 3+q (§))“ n S a n B( fe 2 + a (g ))2 ((fo, ho), cr) such that 

l|9/^r(/*^*)-3/^(/o,/io)|| jC ( A2+a(s)ifcl+0(s) ) < 1/2 ki. 

The perturbation result [2, Theorem 1.1.3.1 (i)] implies 

-df^Ufo, ho) € U ( h 2+a (S ), h 1+a ( S), 2k u Sr), 

which yields the desired claim. □ 


Remark 5.3. We will see in the proof of Theorem 5.7 that, in contrast to K \, the constants ui\ in 
Lemma 5.2 appear to depend on the || • \\o +a -norm of (/*,/i*). Whence, for (/*,/i*) close to (fo,ho), 
these constants may become large. 


We are now left to establish the assumptions of Lemma 5.2 for some sufficiently small a. There¬ 
fore, we pick an arbitrary (/*,/i*) € (/i 3+ “(§)) 2 n S a and introduce a parameter r G [0,1] which will 
enable us to continuously transform the leading order part := df&f (/*, /i*) of <9j < hi(0, (/*, h*)) 

into a (negative) Dirichlet-Neumann map. This will allow us to use the continuation method and 
prove, by relaying on the properties of this Dirichlet-Neumann map, that large positive real numbers 
belong to the resolvent set of df&i(f*, /i*). We emphasize that this construction uses to a large 
extent the additional regularity of (/*,h*), as the mappings (v+,vZ_) introduced in (3.15) possess 
additional regularity close to the boundary To in this case. More precisely, for each r G [0,1] we 
introduce the operator df§f* T G £(/i 2+ “(§), h 1 +a (E>)) by 


Wr[/] 


rk / 2/' tr 0 d y v*_ 
l-i - ^ f*-d 


tr 0 d x vl)f - B(f*)wl T [f], f G h 2+a (§), 


(5.7) 
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with (v* + ,v *_) being defined in (3.15) and (w + T [/], wf T [/]) denoting the solution to 


^(/*AK r [/] = 

^- T [f] - 

= -TdflF{f*,h*)[f]v* + +Td f B*(f*)[f]v*_ 

r , r . . . . , tro tro 

w+ r [/] - [/] = - g(p- - p+) + (1 - t) 


Wlrlf] = 0 

wl T [f] = 0 


h* f * f* d 


in 

n + , 

in 

n_, 

on 

To, 

on 

To, 

on 

r ls 

on 

r_i. 


(5.8) 


As (/*,/>*) € (h 3 +«(§)) 2 ns a , a density argument together with [1, Theorem 9.3] implies that the 
transformed potentials (v+, v *_) satisfy 


tr Q d y v - _ tr Q d y v *+ h 2 +a 
f*~d h, ~ /* 


(5.9) 


and therefore (reT T [f],wf T [f]) is well-defined, cf. Corollary 3.4. For r = 1 we recover the leading 
order part = df<S>J (/*, /i*) of <9/-<3?i(0, (/*, 7,*)), while S/TJq is the above mentioned Dirichlet- 

Neumann map. Our method uses localization techniques based on a suitable partition of unity 
allowing us to keep up the setting of periodic functions. 


Partition of unity. For each p > 3 there exists a family of functions {IT' }l<j<2P+ 1 C C°°(§, [0,1]) 
such that 

(*) supply = U neZ (27rn + /J), where /J := [j - 5/3, j - l/3]vr/2 p ; 

(n) T,f=i = 1 in C(S). 

The center of the interval 7? for 1 < j < 2 P+1 is denoted by x v - := (j — 1)tt/ 2 p . We call the family 
{n^} 1<i<2P +i a p-partition of unity. 

Moreover, let {x*j}i<j< 2 P + 1 C C'°°(§, [0,1]) be a corresponding family of functions such that 
(a) supp x P j = U neZ (27rn + Jj) and /J C Jj; 

(&) Xj = 1 on /J. 

We can achieve that Jj = /J _ 1 U Ij U 7j +1 , where 7 q := I^+i and I^ p+1+1 := If, so that /J and Jj 
have the same center x;^. 

The following remark is a simple exercise. 

Remark 5.4. Given k,p £ N with p > 3 and a € (0,1), t/ie mapping 

/ca max ||II^/|| fc+a 
i<i<2 p+1 J 

defines a norm on h k+a (E>) which is equivalent to the || • ||fc+c«- norm. 

The following perturbation type result stays at the core of our analysis. 

Theorem 5.5. Let a > 0 he such that (fo,ho) € S a and let p > 0 and a! € (0, a) be given. Then, 
given (/*,7*) € (/i 3+ “(§))“ Cl S a , there exist an integer p > 3, a p-partition of unity {IT!?} 1 < ? -< 2 p+i j 











18 


J. ESCHER, B.-V. MATIOC, AND CH. WALKER 


and a constant K 2 = K 2 (p), and for each r £ [0,1] and 1 < j < 2 p+l there are bounded operators 
A T J € C(h 2+a (§) ,h 1+a (§)) such that 


\\rfjdf ,* [/] - Kj [n? /] 11 !+ a < A*||n!?/|| 2+a + at 2 ||/|| 2+q , 


(5.10) 


for all f £ h 2+a (§). The operators A r . j are defined by the formula 


A. 




rk ( 2 /q tr 0 d y v*_ 


f* — d 


- tr 0 <9 x u* 


/'-A 

9 17 


1 + /;■ 


/ 2 


A- ' f* ~ d 


, tr 0 dywli [/] - /o(xp tr 0 d x wfi [/]j, 

(5.11) 


where denotes the solution to the problem 


Ao,j(f*, h *) w +r[f ] = ^+/" mST2 + , 

“ 4 0j(/*) u, -r[/] = r ^-/" inO_, 

#j(/*A)w+r[/] - £j(/*)W-r[/] = r - B / / OB F 0 , 

™+r[/] “W-rE/5 = ~[ A p + (1 ~ t)A a ]/ on r 0 , 

w+rE/j = 0 onTj, 

W-r[/] = 0 on r_i, 


(5.12) 


and 


A + 

B := 


tr 0 cAu* 


h* — /* xA 
k / 2flti' 0 d y vl 
fi-\ f*~ d 


A _ := 


tr 0 


f* — d 


P 

X- 

3 


- tr 0 


k / 2 fjtT Q dyV. 
T+ 


A P := g(p- - p+), A a 
- tr 0 <9 x u)Q 


:= Aj 


d* f* 


A-, 


Moreover, Aoj(f*,h*), A^^f*), Bj(f*,h *), and Bj{f*) are the operators obtained from Au(f*,h*), 
Ao(f*), B(f*,h*), and £?(/*), respectively, when evaluating their coefficients atx p . 

Proof. Let c/ £ (0, a) and p > 0 be fixed. Given p > 3, a p-partition of unity {nj} 1<J<2 p+i, and 
r £ [0,1] we decompose the difference U p df$i* T [f] — A r j[I Tfjf] = T\ + T 2 + T3 by setting 


T\ := 
T 2 := 
T 3 := 


r/c /2/'tr 0 cU;* 


p_ \ /* - d 

k l + f? 


-tr 0 d,v-_)\ja-,fy-U 


rk ( 2/' tr 0 


p_ \ /* - d 


- troa r v*^n^/', 


tr 0 dywff [U p /] - — 1 + II? tr 0 dyW^. [/], 

h- J* ~ a 


P- /* - d 

— /*n? tro dxW_ T [/] - —fi{x p j) tr 0 <9 x u/^[II^/] 

P— J p— J J 


and estimate each of these terms separately. In the following we shall denote constants which are 
independent of p (and, of course, of / £ h 2+a (§>), r £ [0,1], and j £ {1,..., 2 P+1 }) by C. 
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The estimate for T\. Noticing that ||lT’/ / '|| a < || (n 7 ./) // || Q , + TT ||/"||2 and using the fact that 


X^n p - = n P j, we have 


rilll+a <K\\f\\ 1+a + C 


V 


2 /' tr 0 d v v* 


<C 

+ K\\fh- 


pf fltoodyV- 
Xj \ f* ~ d 


V 

x i 


f* — d 
_ fi tr 0 d y v*_ 
f* — d 


- tr 0 d x v*_ 


2 /' tr 0 d y v*_ 
f* — d 


- tr 0 d x v*_ 


mf 


l-\-a 


+ ||x?(tro d x v*_ - «9x^(^,0)|| 0 ||n?/|| 2+ a 


Choosing p sufficiently large and using the uniform continuity of the functions in the square brackets, 
we are led to the estimate 

mi| 1+ c <W3)||n^/|| 2+Q + K\\f\\ 2 . (5.13) 

The estimate for T> and T 3 . The terms 75 and T 3 are estimated in a similar way. However, 
due to the nonlocal character of these expressions, the arguments are more involved than in the 
previous step. It is easy to see that 

ll^Hi+a <C 1 l/*n^ tr 0 dxw’L^f] - f^x?) tr 0 f]\\ 1+Q 

tr 0 <9a; (n^w!5 r [/]) - tr 0 d x w^[H p f] ||i +Q , + K\\f\\ 2+a > 

<C\\x P jK tr 0 d x (n^ r [/]) - X P jfl{x P ) tro d x w*_* [n ?f] ||i +a 

+ C'lK 1 - + K \\fh+a' 

=:T 31 +T 32 + K\\f\\ 2+a ,, 


(5.14) 




where we used once more the fact that y(n( = n. We first note that 

v J J J 

T32 <C \\(1 - x p ) tr 0 4 ^[n^/]|| 1+a < K\\f\\ 2+a , + C||(l - x>_^[nf/]|| 2+Q . 

The pair (u+,u_) := (1 — Xj){ w +t [Hj/]■ W-r [n^/]) is the solution to 

■Ao .■(/*,/i*)«+ = tA+(i - x?)(n P J)" - (x?)"^+r[n?/] 


-2(xp , 4'W+r[n^/] + 


2 /: 


r v 
i 


a, +TL K _ f 

A o,j(f*) u - = rA-(i-$)(!% 


ix^ydyW^in^f] in n+, 


-2(^yd x wii\mf] + 


2 /: 


f*-d 


Xx^YdywlillPf} in Q_, 


Bj(f*,h*)u+ - = tB( 1 - xp( n j/)' + k V+ (X^)7*(^j)tr 0 w+y[n^/] 

-kiiZ 1 {x P )'fi( xP ) tr 0 w*i [n?/]) 
u+-u_ = — [A p + (1 — t)A^]( 1 — xpn^/ 
u+ = 0 
u- = 0 

and, because of (1 — x?)n? = 0 , we end up with 


on r 0 , 

on r 0 , 
on Ti, 
on r_i, 


T 32 < K\\f\\ 2+a , + C\\u-\\%- a < K\\f\\ 2+a ,. 


( 5 . 15 ) 
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The remaining term T 31 can be estimated as 

^31 <C\\x P j{fi ~ fKtf)) tr 0 a x u;^[n^/]]|i +C( + C\\x p tr 0 d x (U p w%[f] - w^[U p f}) ||i +Q 
< (m/6 ) 11 /112+a + K 11 /11 2 +a' + C\\ X P tr 0 d x (U p w n _ T [f] - [U P f])\\ 1+a (5.16) 

if p is sufficiently large. We are left to estimate the last term in (5.16). We note that 

\\x p tT 0 d x (n p w 7 L T [f} -w n ^[u p f])\\ 1+Q <K\\f\\ 2+a , + \\x p {wy T [n p f} -n^ r [/])|| 2 + Q . (5.17) 

We now infer from the fact that the pair 


(z + ,Z-) := X?((^+r[nf/],W-r[n^/]) -iPAwlMwlrlf])) 


solves the diffraction problem 


•A%{f*,h*)z + 


AUf*) z - 




B(f*,h*)z + - 


Z + — Z- 

Z+ 

z_ 


rA + x P (n p f)" - Lg^g+n^/" 

( X p )V4[U p f] - iIi p )"wl T [f] 

^{x P )'d x w^[U p f\ - 2 (U p j )'d x wl T [f} 

hf^(n p )'d y wi T [f] - fef^| xP (xp'^+r[nj/] 

rA-x^f)" - 

(4)V4\n J/] - (nJ)"«£ T [/] 

2(xP'Ato”[n"/] - 2(n 

-ML(W j ya,wi T {f] - jh n?/l 


f*A) - 6 i (/*,/l*))[x^+r[n^/]] 

3(/*) - ^(/*))[x>^[n^/]] + rSx^/y 

- rI^9 / B’ r (/*)[/]«i 

^((Xp'/^^ptro w^ [U p f] - (np7'tr 0 u;^ T [/]) 
/il 1 ((xp 7 * 7 pt r 0 W-r[ n i/] - ( n p 7 *tr 0 wlrlf]) 
r)v p ( TtIodyV + _ Tt *o d v v - _ a )n p f 


0 

0 


in 0 _)_, 


in fl_, 


on T 0 , 

on T 0 , 

on Ti, 
on T_i 
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and the Schauder estimate (3.6) that 

K(/*, M - A^.(/*, K)) [x p jW*4[np}] IU + 

+ \\{AUf*) - ^(/*)) [*j«e?p;/]] iu + 

+ ||(B(/*,M - Bitt-A)) [#i#P?/]] Ili+a 


x J 


,/_tro_5^4 _ 
j Wl,- /* 


A, 


|n */"' 


P / tro 

X H f,-d ~ 


A _ 


|n?/"ll 


+ II (£(/*) - #?(/*)) [X^_r Pj/JJ II l+o 

+ llrBx^f) 1 + TU p j d f B^(f,,K)[f}v* + -TlF j d f B*(f.)[f]v*_\\ 1+a 


+ 


P 


t ti'o d y v+ t tro dy v - 

f* — d 


Aa)|| ||I^/|| 2+ a 


h* f* 

+ K ll/lb+o'- 

Using the same arguments as when estimating T\, we find for p sufficiently large 

11 x P j (W-Tl n j/] - A P jwl T [/]) 11 2 ~ a = \\z-\g~ a < (p/6)\\nPf\\ 2+Q + K\\f\\ 2+a , 

which yields, together with (5.15)-(5.17), that 

m\\ 1 +a <(p/m^fh +a +K\\f\\ 2+a , 

Similar arguments show that 

ll^lk+o <W3)||n^/|| 2+a + K\\f\\ 2+a >. 

Gathering (5.13), (5.18), and (5.19), we have established the desired estimate (5.10). 


(5.18) 

(5.19) 
□ 


Fourier analysis: The symbol of A rj . In this step we use Fourier analysis arguments and ODE 
techniques to represent the operators A T] introduced in (5.11) as Fourier multipliers. Subsequently 
we will use a Marcinkiewicz type Fourier multiplier theorem to prove that all these Fourier multi¬ 
pliers Arj are generators of strongly continuous and analytic semigroups, cf. Lemma 5.6, provided 
the more dense fluid lies below the less dense one and under certain conditions on the initial data. 


We fix now r £ [0,1] and p, j £ N with p > 3 and 1 < j < 2 P+1 arbitrary. Given / £ h 2+c 
consider its Fourier expansion 

/ = E ^e imx 


we 


ra£ Z 


and look for the corresponding expansion of A rJ [/]. In view of (5.11), we needed to determine the 
expansions for tro V {w+l [/], [/]). Hence, we let 


(w+t[/],™-t[/D = X! ( A 

mGZ 


A 

m ’ J 


)fn 


(5.20) 


with functions (A+,A m ) = (A+, A m )(y) still to be determined. Recalling (5.3) and the formulae 
in Appendix A, we use the summation convention and write 

Aoj{f*, K) = c+fipu -4qj(/*) = c~ t d p i, Bj{f *, K) = /3 Z + tr 0 d h Bj{f*) = Pf tr 0 d u 

and we set 


a± ■= cf 2 /cf 2 , b± := 1 /c 22 , D± := s/b± - (a ±) 2 
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noticing that b± — (a±) 2 > 0. From (5.20) and (5.12) we deduce that for fixed m € Z, the pair 
(^4 +,t 4“) solves the following problem 


(A.+ )" + i2ma + (^4+) / — b + m 2 A+ l = —rb+A+m 2 
(^4“)" + i2ma-(A^y — b-m 2 A^ n = —rb-A^m 2 
imffl A+ (0) - /?f i”(0)) + p% (A+ )'(0) - {A~ )'(0) = irmB, 

A+(0 )-A~(0 ) = -[A p + (1-t)A a ], 
A+(l) = o, 


with constants A±, B , A p , defined in Theorem 5.5. We now set 


in 0 < y < 1, 
in — 1 < y < 0, 


(5.21) 


cos±(y) := cos (a±my), 
cosh±(y) := cosh (D±my), 


sin ±{y) := sin(a±my), 
sinh±(y) := sinh(.D±my). 


With this notation, it is easy to verify that the general solutions to the first two equations of (5.21) 
are given by the formula 


At ■■=ut + ivt 


(5.22) 


with real parts 


a± 


u m {y) =£i ( cos±(y) cosh±(y) + — sin± 


+ 


+ 


JL 

D±m 

t 


cos±(y) sinh-t (y) 


^(sin±(y)cosh ± (y) - cos±(y) sinh±(y)) + sin±(y) sinh±(y) + tA± 


and imaginary parts 

v m(v) =€i ( - sin ± (y) cosh±(y) + cos±(y) sinh ± 
+ ^(cos ± (y) cosh±(y) 


D±m 


£ 


D±m 


sin±(y) sinh-t (y) 


1 

+ —sm ± 


+ 


JL 

D±m 


cos ±(y) sinh ± (y). 


The real constants {gt : 1 < i < 4} are to be determined such that the last four equations of 
(5.21) are also satisfied by (A+,yl“). It can be shown by explicit, but tedious computations that 
such constants can be uniquely determined to obtain (.A+,.A“). Since 

tr 0 dywti[f] = ^ (£2 + *£4 )fme imx , tr 0 d x wt 3 T [f] = ^ im(£ 1 + tA_ + i^ 3 )f m e imx , 


mEZ 




and using the definition (5.11) together with the explicit expressions for : 1 < i < 4}, it follows 
that Arj is the Fourier multiplier 


A T ,j[^2fme imx ] = J> mfme ir 
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with symbol (Am)meZ defined by 


Re Am .— 


/tanh(D + m) tanh(D_m) \ - 1 


Pt D + m 


+ 


j5 2 D_m 


A„ + + 


tA_ cos(D_m) rA + cos(D + m) 


cosh (D_m) 


cosh(D + m) 


Im A m :=— 
h- 

+ T 


rk /a-d y v*L(x p -,0) 


+ D 2 _ 


+ d x vf(x p ,0 )jm 


/tanh (D+m) tanh(ZD_m 


P2 D + m 


^ v - 1 'A_sin(.D_m) A_|_ sin(L)+m) 

j3 2 D-m / L cosh(D„m) cosh(L> + m) . 


— T 


/tanh(D + m) tanh(Z)_m) \ - 1 t&nh(D+m) 


Pi D + m 


+ 


j3 2 D_m 


PiD ^ 


■ [Pt A+ — Pi A- — B]. 


(5.23) 


(5.24) 


We refrain from presenting here the detailed computations that are used to determine the constants 
: 1 < i < 4}, (5.23), and (5.24), as they are quite long and only the outcome, that is, the 
explicit formula for the symbol (A m )m£Z, is of importance for the further analysis. 

We are now in the position to prove that the operators A r j are generators of strongly continuous 
and analytic semigroups. To this end we will use a Marcinkiewicz type Fourier multiplier theorem 
[29, Theorem 2.1], which generalizes a result from [4] (see also [19] for a similar result) and states 
that a Fourier multiplier 

£ fme imx ^ £ A m f m e imx 

mEZ mG Z 

belongs to £(C' r+ "(§), C s+a (§)), r, s £ N and a £ (0,1), provided that 


si := sup \m\ s r |A m | < oo 
m€Z\{0} 


and s 2 := sup |m| s r+1 |A m +i — A m | < oo. 

m£Z\{ 0} 


Additionally, the norm of the Fourier multiplier as a bounded linear operator from C r+ "(§) to 
C' s+a (S) can be bounded in terms of the constants s\ and s 2 alone. Bearing this result in mind we 
can derive the following lemma. 


Lemma 5.6. Let cr > 0 be such that ( fo,ho ) € S a . Then, there exist constants > 1 and 
u>\ > 0 depending only on a , such that for any (/*, /i*) £ (/j 3+ “(§)) n S a , any p-partition of unity 
{Ii p } l<j<2P +i with p > 3 and any t £ [0,1], the operators A r j, 1 < j < 2 P+1 , defined by (5.11) 
satisfy 

A - A Tj £ Isom(/i 2+Q (§), h 1+a (8)), (5.25) 

KilKA - A tJ )[/]||i +q > |A| • ||/|| 1+Q + ||/||2 + a (5.26) 


for all f £ /i 2+ "(§) and A £ C with Re A > u )\. 
Proof. We write 


R® A m — p m n m , 
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where 


I'm ■= T 


Hm ■ — 


/tanh(D_|_m) tanh(D_m) \- 1 r A_ cos(IA_m) A + cos(D + m) 

V P 2 D + m P^D-m ) V cosh(D_m) cosh (D+m) J’ 
/tanh(D_|_m) tanh(D_m )\ _1 


+ 


(A p + A a) ■ 


V @2 D+m /3 2 D^m 
Recalling the definition of S a , there exists a constant Co = Co(cr) > 1 such that 

C 0 1 < A p + A^ < Co- 

In fact, it is not difficult to see that we can choose Co large enough to guarantee that 

Cq 1 < -pA < C 0 , m€Z\{0},r€[0,l], 

\m\ 


(5.27) 


sup sup ( \v m \ + 
m€Z\{0} re[0,l] 


Re A m | + I Im A,, 


\m\ 


<C 0 , 


(5.28) 


sup sup (| Re A m+ i - Re A m | + | Im X m+1 - Im A m |) < C 0 . 

m€Z\{0} rG[0,l] V ' 


Let now A = Ci + *£2 € C be given such that Re A = £1 > tui := 2Co- Since Ci — Re X m > Co, the 
(formal) inverse (A — A r j)~ 1 of A — A r j is the Fourier multiplier 


/m 

m£Z 


e imx i-A 


^ Amfn 


with symbol 


» _ / \ 1 _ Ci Re A m *(C 2 Im A m ) 

: — V 7 ' ^ m) — 


(Cl - Re A m ) 2 + (C2 - Im A m ) 2 ■ 

Step 1 . We first prove that (A — A rj -) _1 belongs to £(C 1 + “( 8 ), C 2+ “(§)), the norm being inde¬ 
pendent of p > 3, j G {1,..., 2 p+1 }, r € [0,1], and A € C with Re A > More precisely, we show 
that 1 


sup sup sup sup (|mA m | + |m| 2 |A m+ i — A m |) < 00 . 
Cl>wi x p &> re[o,i] mez* 


(5.29) 


Note that (5.28) implies 
IA-o| < Co, 
|m| 2 |A m | 2 = 


m 


< 


m 


< 


C 2 m 2 


<cl 


(5.30) 


(Ci - Re A m ) 2 + (C2 - Im Xm ) 2 " (Ci - Re A m ) 2 m 2 + C 0 4 ~ 0 ' 
Setting R m '■= (1 — ReA m and I m ■= (2 — Im A m , we have R m > 0 and 

|A m _|_i — A m | < T\ + T2 + T3 + T4, 


lr The operator A t j depends on the partition {n^} 1 < 3 < 2 p+i through the middle point x p of the interval 7? only, 
and therefore A t j makes sense when replacing x p by any 16 S. In view of this fact, when taking the supremum over 
x p € § in (5.29) we consider a larger set than when taking the supremum over all middle point x p and all p-partitions. 
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where 


T\ : = 
T 3 : = 


Rm Rm -\-11 Rm .+1 Rn 

(Rl + Il)(Rl^ +H 


Im I'] 


m 1 m -\-11 l-^rn+1 


C 77l+1 
Im.- 1-1 In 


771+1 ) 


T 2 : = 
T a : = 


\R-rnXm+\ R 


m+i^-rn 


(Rl + ll)(Rl +1+ H 


m+ 1 > 


I T TD 2 _ T p2 


( J R 2 , + / 2 J( J R 2 m + 

-^m+1 )’ {Rm + -Im)(+m+l + l‘rri J r 1) 

The estimates for m 2 T 3 and uXTa are similar to those for m 2 Ti and m 2 T 2 , and therefore we shall 
present only those for the latter. Recalling (5.28) and (5.30), we get 


and 

m 2 T2 <m 2 


I R 


m | \ \ -^7n+l | 


m 2 Ti < m 2 \A r 


+ \Im\)\Im+l R 


l-^-rra+11 

I + l^mll-^m+1 ~ R 


Rm +1 — Rm | ^ 2Cq 


(R 2 m + Il)(R 2 +I 2 ) 


A Cqtti | A m _|_i | (| A m | + 2 |A m+ i|) 


<10 cl 


Proceeding similarly with m 2 T 3 and m 2 T 4 , we arrive at (5.29). In view of [29, Theorem 2.1] we 
additionally know that (A — A r j ) _1 € £(C n+ “(§), C n+1+ “(§)) for all n € N, and since the closure 
of C n+ 1 +Q (S) in C n+Q (S) is exactly h n+a { S), a density argument leads us to the desired property 
(5.25). Moreover, our arguments show that that there exists a constant K\ depending only on Cq 
such that 


ki||(A — A tj )[/]||i +q > ||/|| 2 +a for all / € h 2+a (§) and ReA > wi. (5.31) 

Step 2. We are left to prove that we can choose large enough to guarantee that 
«h||(A - A Tj )[/]||i +a > |A| • ||/||i +a for all / € h 2 +a (S) and ReA > 

For this it suffices to show that 


sup sup sup sup |A| (|A m | + |m| • |A m+ i - A m |) < oo. (5.32) 

Cl >^l x?£§ re [0,1] mez* 


Using (5.28), we first note that 


Cl 2 


< 


Cl 2 


< 


Cl 2 


(Cl - Re A m ) 2 + (C 2 - Im A m ) 2 “ (Ci - ReA m ) 2 “ (Ci/2) 2 
Additionally, (5.28) and (5.30) lead us to 

Cl 


_^ 2(C 2 ~ Im A m ) 2 + 2 (Im A m ) 2 < 

(Ci - Re A m ) 2 + (C 2 - Im A m ) 2 “ (Ci - Re A m ) 2 + (C 2 - Im A m ) 2 “ 

<2(1 + Cq), 


1 + 


= 4. 


Im A m | 2 , 2 1 2 \ 

(m |A m | ) 


m- 


and we thus have established that 

sup sup sup sup |A| • |A m | < 2(2 + Cq). (5.33) 

Ci>^i a;^eSre[0,l] m£Z* 


To hnish the proof, we recall that |A m +i — A m | < Ti + T 2 + T 3 + T 4 . Combining (5.28), (5.30), and 
(5.33), we find that 

|Am|7i < |Am| • |A m | • |A m+ i| • \R m+ i - R m \ < 46+2 + C$) 
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and 

| \m\T 2 <Co|Am| • |A m _|_i| (|A m | + 2|A m _(_i|) < IOCq (2 + Cg). 

Proceeding similarly with the remaining terms |Am|T3 and \Xm\T 4 , we obtain (5.32). Together with 
(5.31) we deduce (5.26), and the proof is completed. □ 

With these preliminary results, we can now prove that the assumptions of Lemma 5.2 hold and 
establish in this way (3.16). 

Theorem 5.7. Let (/o, ho)) € V and 60 € h 2+Q (§) be given such that the first inequality in (2.7) is 
satisfied. Then 

0, (/o, ho)) € n{h 2 + °(S), h 1+ “(§)). 

Proof. Recalling (5.4) (with (/*, h*) = (/o, ho)) and the perturbation result [2, Theorem 1.1.3.1 (ii)], 
we are left to show that — <9y<3?^(/o, ho) € "R(h 2+ "(§), h 1+ "(§)). To this end, it suffices to verify the 
assumptions of Lemma 5.2 for some a > 0. 

Let a > 0 be such that (/o, ho) £ S a , let and be the constants found in Lemma 5.6, and pick 
of £ (0, a). Given (/*, h*) £ (h 3+ "(S)V n S a , Theorem 5.5 implies the existence of a p-partition of 
unity {n^} 1<i<2P +i and of a constant iL 2 such that 

l|np/$p* r [/] - A T;j [n^/]|| 1+0! < ^-||n^/|| 2+a + iL 2 ||/|| 2+Cf / (5.34) 

for all / € h 2+ “(§), r € [0,1], and j £ {1,..., 2 P+1 }. In view of (5.26), we get 

«r||(A - A Tj )[n^/]|| 1+Q > |A| • ||n?/|| 1+a + ||n*/|| 2+a (5.35) 

for all A £ C with Re A > uq, r € [0,1], / € h 2+ "(S), and j £ {1,... ,2 P+1 }. Combining (5.34) and 
(5.35) gives 

«!||d?(a - d f ^)\f]\\i+ a >«il|An y - a tJ [i^/]|| 1+q - Ki \\a rJ [n^/] - n?0/<i>r* [/]|| 1+Q 

>|A| • ||nf/|| 1+Q + unf/lla+c - |||n^/|| 2 +« - *4if 2 ||/|| 2+( y, 

hence 

KlK 2 \\fh + a' + «r||I^(A - df$l* T )[f]\\i +a >|A| • ||n^/|| 1+Q + ^WlSPjfh+a 

for all Re A > uq, f € h 2+a (§), r £ [0,1], and j £ {1,... ,2 P+1 }. Together with Remark 5.4, we 
conclude that 

^AAH/lh+a' + Kill (A - df*ir)[f] llr+a >|A| ■ ||/||l+a + |||/|| 2 +« 

for all Re A > u\, f £ h 2+ "(§), and r £ [0,1]. Using the interpolation property (4.9) and Young’s 
inequality, we find a constant uq > 0 depending on R 2 such that 

4*4 || (A - )[/] || 1+ a >|A| • ll/lll+a + II/II 2 +a (5.36) 

for all A € C with Re A > uq, / £ h 2+a (S), and r £ [0,1]. Additionally, choosing o sufficiently small 
we obtain from the definition of S a that 

hi := max |4 ki, sup {||9/$^)|| 2 + a(S))/ll+ „ (s) ) : (/*>M e D %2+« (s)) 2((/ 0 , h 0 ), a)}} < oo. 
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Setting r = 1, we conclude that if (cui — = (wi — df$J (/*, h*)) G fsom(/i 2+a (§), h 1+ “(§)), 

then 

-d f $l (U,K)en (h 2+a (S), h 1+a (S), ki , wi), 

where is independent of (/*,h*) € (/i 3+ “(§)) 2 H 5 CT n ((/o, ho), <r). 

Hence, the assumptions of Lemma 5.2 hold true for sufficiently small a if we can show that 
(uj i — G fsom(h 2+a (§), h 1+ “(§)). Using the method of continuity, cf. e.g. [26, Theorem 5.2], 

and (5.36) it follows that (A — is onto for Re A > u> i, provided that (A — is onto. 

But the surjectivity of (A — df&J q) for positive A follows as in the proof of Proposition 5.8 below, 
and we thus have verified the assumptions of Lemma 5.2. □ 


We present now a generation result for the Dirichlet-Neumann operator df&f q in a slightly more 
general context. 


Proposition 5.8. Let (/*,/i*) € V be given and a G h 2+Q (§) be a positive function. The linear 
operator 

A[/] :=£(/>%[/], feh 2 +a ( S), 

with (w 7 L 0 [f],wf 0 [f]) denoting the solution to 


satisfies 




-4o if*,h*)wl 0 [f] 

= 0 

in f 2 + , 

AUf*)^- 0 [f] 

= 0 

in 

B(/^*Ko[/]-fi(/>- 0 [/] 

= 0 

on T 0 , 

“+ 0 [/i - w -oi/i 

= af 

on T 0 , 

*Uf\ 

= 0 

onT i, 

w-oif] 

= 0 

on T_i 


-A G 77(h 2 + “(S),h 1+ "(§)). 


Proof. A short inspection of the proof of Theorem 5.5, Lemma 5.6, and Theorem 5.7 reveals that 
there exist constants k > 1 and oj > 0 such that 


«ll(A-A)[/]|| 1+Q >|A|.||/|| 1+ a + ||/|| 2+ a 

for all A € C with Re A > u> and all / G h 2+a (S). To finish the proof, it suffices to prove that for all 
positive A, the operator (A — A) G £(/i 2+ “(§), h 1+ “(§)) is onto. Let thus A > 0 and F G h 1+a ( S) 
be given, denote by (z+,Z-) G /i 2+ "(fl+) x /i 2+ "(fl_) the solution to the diffraction problem 

Aq (f*,h*)z+ = 0 
•Ao(f*)z- = 0 
B(f*,h*)z+ - B{f*)z- = 0 
-A a~ l (z+— Z-) + B(f*)z- = F 

z + = 0 
z- = 0 

and set / := — a -1 tro(z+ — zJ) G /i 2+ "(§). ft is easy to see that (w+g, W-o)[f] = ( z +i z -)- Therefore, 
(A - A)[/] = A/ + B(U)w%[f] = -A a-\z+ - zJ) + B{f*)z- = F. 


in U + , 
in n_, 
on T 0 , 
on T 0 , 
on Ti, 
on T_i, 


(5.37) 
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In the remaining part of the proof, we establish that for each F E h 1+ “(§) and A > 0, the diffraction 
problem (5.37) possesses a unique solution (z+,z—) € h 2+a (Q+) x h 2 +a (fl_). In fact, due to the 
arguments in the proof of Theorem 3.2, it suffices to prove that the mapping 


( 


(z+,z-) ^ 


\ 


Alo(/*,/i*)^+ \ 

(f*)z- 

B{f*,h*)z+ - B(f*)z- 
Xa 1 tr 0 (2+ — Z-) + B(f*)z- 
tri z + 
tr_i 


{5.38) 


is an isomorphism between C 2 +a (Ll + )xC 2 +a {Ll-) and C' Q (f2 + )xC'"(n_)x (C' 1+a (§)) 2 x ((7 2 + 0 ! (S)) 2 . 
It is easily seen that for A = 0 the operator (5.38) is an isomorphism between these spaces (as the 
third and fourth operators defined by (5.38) lead to decoupled equations). As the mapping 


(z+, Z- ) eA (0,0, 0, —Aa 1 ti(){z + — Z-), 0, 0) 


is compact, we conclude that (5.38) defines a Fredholm operator of index zero. We are left to show 
that (5.38) defines for each A > 0 an operator which is one-to-one. So, let (z+, z~) be a solution to 
(5.37) corresponding to F = 0 and assume that 


max2_ = z_(xo,0) > 0. 

n_ 


If z_ ^ 0, Hopf’s lemma ensures that 

(B(f*)z-)(x 0 ) = A(Ap) _ 1 [z + (x 0 ,0) - z_(x 0 ,0)] > 0, 


hence max^ > max^ z_. On the other hand, if max^ z + = z+(xi,0), then Hopf’s lemma 

implies {B(f*,h*)z+)(x i) < 0, whence {z + — z_)(xi,0) < 0. This contradicts the inequality 
rriaxQ z + > max^ z_, meaning that = 0. But then also z + = 0, and the proof is complete. □ 


6. The second diagonal operator 

In this section we prove that the Frechet derivative ^$ 2 ( 0 , (/o, ho)) is the generator of a strongly 
continuous and analytic semigroup as stated in (3.17) when (/o, ho) and &o are such that the second 
inequality of (2.7) is satisfied. To derive the corresponding Theorem 6.5 we proceed in a similar 
way as in Section 5 and first identify the “leading order part” (/01 ho) of < 9 ^* 1 * 2 ( 0 , (/o,ho))- In 
Lemma 6.1 it is shown, however, that the latter is related to the solution operator of a Dirichlet 
problem which differs from the case considered in Section 5 where the leading order part df&J was 
related to a diffraction problem. The arguments that follow are thus somewhat simpler than those 
in Section 5. 

Lemma 6.1. Let (/*,h*) € V and let h*) € C{h 2+a (E>), h 1+a (§)) denote the operator 

defined by 

/ 2 / 7 ^ r) 77 * \ 

<«(/*, h*)[h] := ( * 1 , y + - tn d x v%)h' - BM*, K)Wl[h\ , h € h 2 + “(S), (6.1) 

ii + \ h* - /* / 
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where (v+,v*L) := (v*/,v*_)(f*,h*,bo) is defined in (3.15) and where W+[7i] denotes the solution to 
the Dirichlet problem, 


Af(f*,h.)W$[h] = 


K - /> 

Wf[h} = 0 
Wm = gp+h 


tr i dyV+ h" in n+, 


on r 0 , 
onT \, 


with 


A\{f*,K) := d xx - 


2 hi 


~9 X y + 


h? + 1 


:d v 


( 6 . 2 ) 


(6.3) 


K-u y (ft *-/*) 2 '" yy - 

Then, given e G (0,1), there exists K 3 = K^e) > 0 such that 
113^2(0, (/*, h*))[h] - &$£(/*, K)[h] ||i +Q < e\\h\\ 2+a + K 3 \\h\\ 1+a for all h G h 2+a (S). (6.4) 


Moreover, d h <f>% € C U (V, £(/i 2+ “(§), /i 1+ "(§))). 

Proof. The regularity assertion follows by using Theorem 3.2. As for (6.4), let e € (0,1) be given 
and a! G (0, a) be fixed. Given 5 G (0,1), we pick a cut-off function x := x& € C'°°([—1,1]) such 
that 0 < x < 1, X = 0 f° r V < 1 — <5, and x = 1 for y > 1 — 5/2. We then obtain from (6.1), (3.12), 
and the Appendix that 

114^2(0, (f*,h.))[h]-d h $Z(f',h.)[h ]|| 1+Q (||^-111-pet + \\Bi{f*,h*)(W + [h\ — W+[/i])||i+ a ) 

<G(||/ l || 1+Q + ||tr 1 ^(W + [/ l ]-^[/ l ])|| 1+Q ) 

<C{\\h\\ 1+a + || X (W+[h] - Wl\h\) ||?+J, 


where (W + [h],W-[h\) = (dhV+(f*,h*,bo)[h\,dhV*_(f*,h*,bo)[h\) is the solution to (3.14) and C is 
independent of 5. We now notice that the pair (tt_|_,tt_) := (Wf[h],0) — (W+[/i], W_[A]) solves 
according to (6.2), (3.14), and the formulas for A(f*,hfi) and df 1 A(f i ,, h*) from the Appendix a 
diffraction problem of the form 


Ai(f*,h*)u+ = Qq h + afh' + (yd y v+ — tri dyvffia/f h n 




+(y - 

l)(^o 9 xy 

W+[h] + bfdyyW + [h}) + btdyW + [h} 

in 

n + , 

•4(/* 

}u- -- 

= 0 



in 

n_, 

£>(/* 

h- -- 

= d h B{U 

,K)[h]v* + 

+ B(f*,K)W$[h] 

on 

r 0 , 

u + - 

U- - 

= 0 



on 

To, 


u + = 

= 0 



on 

Ti, 


It- = 

= 0 



on 

r_i 


with af,bf G h a (Q + ), 0 < i < 2. Therefore, (u+jvff) := x( u +i u - ) solves 

•4T(/*)M«+ = xAlif*, h*)u+ - j^-x'd x u + + (2x'dyU+ + x"u+ ) in fi+, 

A(f*)u *: = 0 in fl_, 

&{f*,h*)u\ - B(f*)u* = 0 on T 0 , 

u\ — wf = 0 on To, 

u^/_ — 0 on Ti, 

wf_ = 0 on T_i 
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Recalling the Schauder estimate (3.6), we obtain that 

KllSa <C\\ X AUU,K)u + \\^ +C(5)\\u + \\^ a , 

<C{\\x(y~ l)(b£d xy W + [h] +btdyyW + [h])\\%+ + \\x(yd y v* + - tr x d y v* + )a$ti'\\%+) 

+ C(5)\\h\\ 2+a , 

<C(\\x(y~ l)|lo + + II X(yd y v+ ~ tri9 y ?;+)||o + )||h||2+a + C(5)\\h\\ 2+ a', 

and choosing 5 > 0 such that \\x{y ~ l)|lo + + \\x(y^y v + ~ tr i ^y v +)llo + < e/2 C, interpolation 
properties of small Holder spaces lead us to the desired estimate (6.4). □ 


As in the previous section we introduce for a > 0 the set R a consisting of those (/*,h*) €E V 
satisfying the inequalities 

(1) a < min{/* - d, h* - /*}, \\f*h + \\ h *h<cr~ 1 , 

( 2 ) + 

(3) || tri d y v* + ||o < cr -1 , 


where (vf,vf) is defined in (3.15). Since the functions (/o,ho) and &o are chosen such that the 
second inequality in (2.7) is satisfied, we may choose a such that (/q, ho) € R a . 


Lemma 6.2. Let a > 0 be such that (/q, ho) € R a . Assume that there exists a constant k 2 := k 2 (a) 
and for each (/*,h*) € (h 3+Q (§)) n R a PI M( h 2+ a (§^2 ((/o, ho), cr) there exists a further constant 
oj 2 > 0 with the property that 

-d h ^ 2 (U, K)£U (h 2+ “(§), h 1+Q (§), k 2 , u 2 ). 

It then holds 

-d h <mfo, ho) € 77(h 2+Q (§), h 1+ “(§)). 

Proof. The proof is similar to that of Lemma 5.2. □ 


We are thus left to prove that there exists a > 0 such that dh & 2 (/*, h*) is a analytic generator 
for each (/*, h*) € (h 3+ “(§)) (^2+c(§))2 ((/o, ho), cr), with the constant k 2 depending only 

on cr. To this end we use the additional regularity of (/*,h*) € (h 3+ “(§)) 2 n to introduce 
again a parameter r £ [0,1] which enables us to continuously transform the leading order part 
dhQlf* := (9/j<!>!(/*, h*) into a (negative) Dirichlet-Neumann map. More precisely, for each r € [0,1] 
we define the operator dh& 2r € £(h 2+a (S), h 1+ “(§)) by the formula 

/ 2/7^ fri \ 

4^; T [h] :=-( * J + -tr 1 4<)h / -g 1 (/„h,)^ T [h], h € h 2+ “(§), (6.5) 

/i-|- \ /* ' 

with HT r [h] denoting the solution to 


Af(f*,h*)Wf r [h\ 
Wf T [h] 
W+ T [h] 


r tri dyV+ ^n 
h* f* 

0 


hh+ - (1 


tri <9^+- 

h* y* - 


h 


in f2+, 
on T 0 , 
on R. 


( 6 . 6 ) 
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For r = 1 we see that dh& 2 *i = dh ^ 2 (/*> h*), while for r = 0 we obtain a Dirichlet-Neumann 
operator. 

We now prove the following perturbation result. 


Theorem 6.3. Let a > 0 be such that (fo,ho) E R a and let p > 0 and a' E (0, a) be given. Then, 
given E (/i 3+ "(§))^ C\R a , there exist an integer p > 3, a p-partition of unity {nf’} 1<J<2 p+i, 

and a constant K 4 = K^{p), and for each t E [0,1] and 1 < j < 2 p+l there are bounded operators 
B T j E T(h 2+a (§), h 1+ "(§)) such that 

WU^dh^lh] - B rJ [n^]|| 1+a < n\\n%h\\ 2 +a + K4h\\2+a> (6.7) 

for all h E h 2+a (E>). The operators B r j are defined by the formula 


™ n.i rk (2h' Jf ti l d y v% * 

B rj[h\ := -( —7 - - -tri d x v A 


T+ 


h* f* 


. k ( 1 + h ' 2 
h - 


f* 


where Wff [h] denotes the solution to the problem 

t tri d y v* 


tri dyWl’*[h]-h*(x?) tri d x Wlf[h]j , 

( 6 . 8 ) 


Al^KWZfih] = - 1 -^_ 
Wli[h] = 0 * 


wH\h] = g P+ -(l-r) 


tri d v v* A 


h* f* 


in n_|_, 

on r 0 , 

h on Tj, 


(6.9) 


with Afj(f*, K) being the operator obtained from AJ(f*,h*) when evaluating its coefficients at x p . 

Proof. Let p > 0 and a' E (0, a) be fixed. Given p > 3, a p-partition of unity {n^} 1<J<2 p+i, and 
an associated family {x J j}i<j<2P + 1 ( see Section 5) we decompose 

U p d h ^ 2 * T [h] - MrjiUPh] = 5i + S 2 + S 3 


with 


Si ■= (1 + T)k9P+ (n p Kh> - K(x p )(u p h)'), 


s 2 ■■= 


S, : = 


T+ 

(1 + r)k 
T+ 
k 1 + K 2 


K tri d v v* A 


h * /* 


(W- 


x p 3 
3 


KtXldyVf ■ 

K - /* J . 


/^+ f* 


tudyWfg^h]- 


k 1 + h ' 2 
p+ h* — /* 3 


nftn d y w«[h]. 


Still, we use C for constants which are independent of p, constants depending on p being denoted 
by K. Recalling that y^TL = II p , we estimate Si as 

IlSilli+a <c\\n p Kti - ^(^)(n^)'||i +a < c\\j\ p h' x p {K - K(a*))\\ 1+a + K\\h\\ 1+a 
<C\b$(K - K(x P j))\\o W^hh+a + K\\hh, 

and similarly 


WS 2 W 


1+0 Y 


<C 


p 

Xj 


2/i( tri d v v*. 


h* f * 


2 h ‘(tri d y v* A 
h * /* 


o l|n^,|| 2+Q + iqh|| 2 . 
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Choosing p sufficiently large, we find that 


ll-Srllr+a + ||S 2 ||l + a <f 11^112+a + K\\h\\ 1+a . 


( 6 . 10 ) 


We are left to estimate S 3 . For this we note that 
i + K 2 


1 +a <C 

<c 

+ 


K - f* x? 


1 + h' 2 

trr d y Wli[Tl p .h\ - trr d y W” T [h] 


(i-4) 


1 + h' 2 


3 K - f* 


h* - /* 3 

MdyWliilPh} 


1+a 


X 


1 + K 2 


3 K - f* 


1+a 

1 + h' 2 


trr dyWli^h] - X P jT —j- trr d y {l^Wl T [h]) 

>1* J* 


1+a 


= :S 31 + S 32 , 


where 


ll^srllr+a <C|| tn d y ((l - ^Wli^h]) | i+Q < C||(l -4)W+ J [n^]||^. 


As z+ : = (1 — Xj) W+r [+h] solves according to (6.9) and (6.3) the problem 

•Aij(f*>h*)z+ = 


t tri d y v* 


K - f* x p 


-2{ x ]yd x wli[n]h] + 


2h' 


h* ,/l 


z + = 0 


z+ = 


9P+ - (1 - t) 


ti'i dyV*> 




(! - x?)n )h 


^ySyWli^h] in n + , 

L j 

on r 0 , 
on Ti, 


we infer from (1 — +) 1+ = 0 and Schauder estimates for elliptic Dirichlet problems that 


IlSnllr+a <K\\h\\ 


2+a' 


( 6 . 11 ) 


Finally, 

||5 3 2|| 1+ a <C 
and for p large enough we have 


p Z' 1 + 4 2 

1 + 4 2 \ 

x H/ 4 -.A 

A h* — /* / 


m\2+ a +c\\x p Awiimh\ - mwiM) fc, 


X 


1 + A 


/2 


3 V /t* — /* +) 


1 + A 


/ 2 


^ h* — /* 


< -+. 

0 - 4C 


( 6 . 12 ) 
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On the other hand, u+ := Xj(W^[Ujh] — solves the Dirichlet problem 

AlU*,K)u+ = K) - K)) [tfWli[mh]] 


and therefore 


+ - 


r tri d v v*. 


h* f* 


xmw - 




h* f* 


H^) H Kr\^h] - (n tyw^ih] 
+2tf j ya x w!g\n? i h] - 2{myd x wi T [h] 


2 h' 


2 h' 

(X P j ydyWl^h} + —^- 




u + = 0 


ll«+ll£« <C 


{^)'d y Wl T [h] in0 + , 

on r 0 , 
on Ti, 


+ 


+ 


y*i (/., K) - Aljtf., K)) [tfwgvqh]] 

jf tri d y v* + tri d y v\ ^ 


X 


X 


j V K - /* 

p f tri d y v* 
1 V K- /* 


h* /* 

tri 


/i* /* 


lin^n 


2 + 0 ! 


+ K ll^lb+o'- 

For p sufficiently large, we obtain together with (6.12) that 

-A* http 


11 <S , 32 11 l+Q! < 


. M. 

4 J 


n + ^ ^ http 


mh\\ 2+ a +K\\h\\ 2+a '. 


2 +a + C'||'U+|| 2 + Q: < -||ilj+i||2+o 

Gathering (6.10), (6.11), and (6.13), the desired estimate (6.7) follows. 


(6.13) 

□ 


Fourier analysis: The symbol of B rj -. Let r £ [0, 1], p, j € N with p > 3 and 1 < j < 2 P+1 be 
arbitrary. Given h £ h 2+a (§>), we consider its Fourier expansion 

h=Y J h m e imx 

mGZ 

and look for the corresponding expansion of M T j[h\. Recalling the definition of M t j from (6.8), we 
first determine the expansions for tri V1F+^ [h]. Let 

Wli[h\ = Y, Bmh m e imx (6.14) 

mGZ 

with functions B m = B m (y ) to be determined. Recalling (6.3), we let (using the summation 
convention) 

A-ij (/*,h*) = Cpid p i 


a := C 12 /C 22 , 6:=1 /c 2 2 , D := \/b - a 2 


and set 
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noticing that b — a 2 > 0. Since W^_f[h] is the solution to (6.9), it follows from (6.14) that for each 
m G Z, the function B m solves the problem 

(B m )" + 2mia(B m )' — bm 2 B m = —rbVm 2 

Bm{ 0) = o, 

B m ( i) = gp+ - (i - t)V, 


in 0 < y < 1 , 


where 


V : = 


tri d y v + 


x p 

J 


(6.15) 

(6.16) 
(6.17) 


h* /* 

The general solution of the first equation of (6.15) is given by the formula 

m •— T iVm 

with real part 

u m (y) =Ci ( cos (amy) cosh (Dmy) + — sin (amy) sin (Dmy)^J + -jj— cos (amy) sin (Dmy) 

+ £3 ( sin(amy) cosh(Hmy) — — cos (amy) sin (Dmy)^j + sin (amy) sin (Dmy) + tV 


and imaginary part 

v m (y ) =Ci ( — sin(amy) cosh {Dmy) + — cos (amy) sin (Dmy)^j — jy— sin (amy) sin (Dmy) 
+ C 3 ( cos (amy) cosh (Dmy) + — sin(amy) sin (Dmy)^j 


Dm 

C 2 _ 

~)m 

C 4 


+ 


Dm 


cos (amy) sin (Dmy). 


The real constants : 1 < i < 4} are to be determined such that the last two equations of (6.15) 
are also satisfied by B m . Explicit computations (much simpler compared to those in Section 5) show 
that such constants can be uniquely determined to obtain B m . Furthermore, using the expressions 
for Q it follows again by explicit computations that the operator B T)/ - from ( 6 . 8 ) is a Fourier multiplier 
with symbol (<^ m ) m , the real part being given by 

Re ip m •— Pm (6.18) 

where 


l^m • — 


k 9P+-V 


fi + tanh (Dm)/m 
The imaginary part is given by 


and 


• — 


Im tp m = 


rk 

P+ 


q((2 - t)V - gp+) 
D 


m + 


rkV cos (am) 
sinh (Dm)/m 

V sin (am) 


sinh (Dm)/m. 


(6.19) 


The representation of B r .j allows us now to prove the following generation result. 


Lemma 6.4. Let a > 0 be such that (/q, ho) G R a . Then, there exist constants k 2 > 1 and w 2 > 0, 
depending only a such that for all (/*, h*) G (h 3+a (§)) nR a , any p-partition of unity {IIj} 1 <y < 2 P+i, 
p G N with p > 3, and any r G [0,1], the operators B T j, 1 < j < 2 P+1 , defined by (6.8) satisfy 

A - M tJ G Isom(h 2+a (§), h 1+ "(§)), (6.20) 

« 2 1| (A — B TJ -)[h]||i +a > |A| • ||h||i +Q , + ||h|| 2 +a (6.21) 

for all h G h 2+ “(§) and A G C with Re A > cu 2 . 
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Proof. Recalling (6.16) and the definition of R a , we find a constant Cq > 0 depending on a such 
that 


Co " 1 < gp+ — V < Co, 


and such that the relations (5.28) are satisfied when replacing (A m ) by (tp m )- With this observation, 
the desired result follows along the lines of the proof of Lemma 5.6. □ 


We are now in a position to establish (3.17). 

Theorem 6.5. Let (fo,ho)) € V and bo £ /i 2+ “(§) be given such that the second inequality in (2.7) 
is satisfied. Then 

-d h * 2 (0, (/ 0 , h 0 )) £ U(h 2+a (S), h 1+a (S)). 

Proof. Recalling (6.4) (with (/*,/i*) = {fo,ho)) and [2, Theorem 1.1.3.1 (ii)], we are left to show 
that —dh^ 2 {fo,ho) £ %(/i 2+ “(§), /i 1+ “(§)). For this, we prove that the assumptions of Lemma 6.2 
are satisfied provided that a is sufficiently small. 

Let a > 0 be such that (fo,ho) £ R a , let K 2 and LO 2 be the constants found in Lemma 6.4, and 
pick a! £ (0, ct). Given (/*,/i*) £ (/i 3+ "(§))^ n R a and using the same arguments as in the proofs 
of Theorem 5.7, Theorem 6.3, and Lemma 6.4, it follows that there exists a constant u >2 > 1 such 
that 


^ 2 \\{\-d h ^2* T )[h]\\ l+a >|A| • \\h\\ 1+a + \\h\\ 2+a 


( 6 . 22 ) 


for A € C with Re A > U2, h £ /i 2+a (§), and r £ [0,1]. As in the proof of Theorem 5.7 we may 
now choose a small to find a constant P 2 > 1 depending only on a such that, for (/*, hfi) belonging 
additionally to the ball B^ 2 + Q re)) 2 ((/o, ho), a), we have 

-d h <S>Uf*,K)eH (h 2+a (§), h 1+a (§), k 2 , w 2 ), 

provided ( 0 J 2 — 9^4? 21 ) = (^2 — 9/i4>J(/*, h*)) is an isomorphism. 

We now establish this last property. Note that, due to (6.22), (£2 — Qh $fi) is an isomorphism if 
( 0 J 2 —52^2,0) i s on to. To prove the latter let A > 0 and H £ h 1+a ( S) be given. We let z € /i 2+ “(fl+) 
be the unique solution to the elliptic boundary value problem 


/ tri d v vl \~ 1 
A [gp+ - h j z + Bi(f*,h*)z = H on Ti. 


Af(f*. hfiz = 0 in 0_|_, 
z = 0 on To, 


(6.23) 


We then set 


K ~ fi 


h ■= gp+ - 


tri d v vl\ -1 


. t r lZ £h 2+a {§). 

ll* J* ' 

Recalling the definitions (6.5) and ( 6 . 6 ), it follows that WT 0 [^] = z and 

(A - d h ^ 2 *o)[h] = A h + Bi(/*, h*)Wlo[h\ = A (gp+ - tfl ^+ )~' tri z + £h(/*, h*)z = H. 


In view of Lemma 6.2 the proof is complete. 


□ 


Finally, we are ready to prove Theorem 2.1. 
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Proof of Theorem 2.1. Because of the equivalence of the problems (2.3) and (3.1) and the reduc¬ 
tion of the latter to (3.8), we are left to investigate (3.8). For the existence and uniqueness result 
we shall employ an abstract result for fully nonlinear problems [30, Theorem 8.4.1], Let us first 
note that (3.11) is implied by Lemma 4.1, Theorems 5.7, 6.5, and [2, Theorem 1.1.6.1 and Remark 
1.1. 6 .2]. Thus, the assumptions of [30, Theorem 8.4.1] are satisfied owing to (3.10), (3.11), and the 
interpolation property (4.9). This proves the existence and uniqueness part. 

Finally, the continuous dependence of the solution on the initial data follows from [30, Theorem 
8.4.4], □ 


7. The Muskat problem with surface tension effects 


In this section we investigate the Muskat problem introduced in Section 2 when allowing for 
surface tension effects in the presence or absence of gravity. More precisely, instead of being con¬ 
tinuous along the interfaces T(/) and T(/i), we assume that the pressure jump along an interface is 
proportional to the curvature of the respective interface, i.e., the pressure obeys the Laplace-Young 
equation 

P-~P+ = -7/«(/) on T(/), 

P+ = ~7 hn(h) on T(h), 

where 7 j [resp. 7 ^] is the surface tension coefficient at the interface T(/) [resp. T(/i)] and where for 
each ( £ C 2 (S) the function 

(" 

/v(U (1 + C ' 2 ) 3 / 2 

is the curvature of the graph [y = ("]. 

Hence, instead of (2.3a) we consider the system 


A u + 

= 0 

in 

n(f,h) 

A u- 

= 0 

in 

«(/), 

d t h 

= — fc/t+ i \/l + h' 2 d v u + 

on 

m, 

u+ 

= gp+h-^ih) 

on 

m, 

U- 

= b 

on 

r d> 

U + — U- 

= g(p+- p~)f + r yfK(f) 

on 

r(/), 

dtf 

= -k/J^y/l + f 2 d v u± 

on 

r(/), 


(7.2) 


supplemented with the initial conditions (2.3b). 

The main result regarding problem (7.2) is following theorem. Its proof is given at the end of 
this section. 


Theorem 7.1. Let g > 0, (/o,ho) £ V (7 (h 4+ “(§)) 2 , and b be given such that (2.4) holds. Then, 
there exist a maximal existence time Tq := To(/o, ho) £ (0, T] and a unique classical Holder solution 2 
(/, h, u_|_,u_) of (7.2) and (2.3b) on [0, To). Additionally, the solutions depend continuously on the 
initial data. 


In order to prove Theorem 7.1, we first recast the problem as a nonlinear and nonlocal evolution 
equation. To this end, we infer from the proof of Theorem 3.2 that for each pair of functions 
(/, h) £ V PI (h 4+ “(§)) 2 and b £ h 2+a (E>) there exists a unique solution 

(v+,v-) := (v+(f,h,b),v-(f,h,b)) £ h 2+a (Q + ) x /i 2+ "(H_) 

2 The notion of classical solution for (7.2) and (2.3b) is the same as in Section 2 with the modification that we 

require additionally to (2.5) that (/, h) G C([0,To), (/i 4+a (§)) 2 ). 
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of the diffraction problem 

A{f, h)v+ 
•A (f)v- 
t B(f,h)v+-B(f)v- 

V+ — V- 

V+ 


0 

in 

0 

in fi_, 

0 

on r 0 , 

9(p+ ~ /»-)/ +7/«(/) 

on r 0 , 

9P+h - 7 hn(h) 

on Ti, 

b 

on r_i 


(7.3) 


the mapping 

[(/,M) ^ (v + (f,h,b),v-(f,h,b))] G C“((vn (h 4+ "(§)) 2 ) x h 2+a (§),h 2+a (n + ) x h 2+a (n_)) 

being real-analytic as a consequence of k € C' a '(h 4_l ~ a (§), h 2+a (S)). Using this observation, we 
deduce that the problem (7.2) and (2.3b) is equivalent to the evolution equation 


= ( f,h )), 

where : [0, T) x (Vfl (h 4+Q (§)) 2 ) C R x (h 4+ "(§)) 2 — > (h 1+ “(§)) 2 is the operator 
defined by 

:= —B(f)v-(f,h,b(t)), 

*2(t,(f,h)) :=-Bi(MK(M,6(i)), 

and where (u_|_,u_) denotes now the solution operator for (7.3). Since we have 

$eC([o,T) x (vn(/i 4 + "(§)) 2 ),(/ i 1 + “(S)) 2 ), 

$(V) G C w (vn (h 4+ “(S)) 2 ,(/i 1+ "(§)) 2 ) for all t G [0,T), 


(7.4) 

($i,$ 2 ) 


(7.5) 


(7.6) 


our next goal is to show that 

- (/o, h 0 )) G ^((h 4+ “(§)) 2 , (h 1+ “(§)) 2 ). (7.7) 

In the remaining part we let (v+,v-) denote the solution to (7.3) determined by the tuple 
(/o,/io,&o) with feo := 5(0). As in the first part, the derivative (/o, ho)) is a matrix 

operator, three of its components being given by (3.12) (with (/*,h*) replaced by (/o,ho)), but 
with 

(w + [f],w~[f]) := (d f v + (fo,ho,bo)[f],d f v-(fo,ho,bo)[f]) 
being the solution to the diffraction problem 

A(f 0 ,h 0 )w + [f] = -dfA(f 0 ,ho)[f]v+ 

A{f 0 )w-[f\ = -d f A(fo)[f]v- 

B(f 0 ,ho)w+[f]-B(fo)w-[f] = -dfB(f 0 ,ho)[f}v+ + dfB(fo)[f]v- 
™+[f] ~ w -lf} = 9(p+ ~ P-)f + 7/3/k(/o)[/] 


w+[f] 

w-if] 


= 0 
= 0 


and with (W + [h\,W-[h\) := (d h v + (f 0 ,ho,bo)[h\,d h v-(f 0 ,h 0 ,b 0 )[h]) solving 


A(f 0 ,ho)W+[h] 

A(f 0 )W-[h] 


= -d h A(f 0 ,h 0 )[h\v + 
= 0 


B(fo,ho)W + [h]-B(f 0 )W-[h] = -d h B(f 0 ,h 0 )[h]v + 


W + [h} - W-[h] 
W + [h] 
W-[h} 


= 0 

= gp+h-^hdh^{ho)[h\ 
= 0 


in 

n + , 

in 

n_, 

on 

r 0 , 

on 

r 0 , 

on 

n, 

on 

r_i 

n + , 


n_, 


r 0 , 


r 0 , 


Ti, 


r_i. 


(7.8) 


(7.9) 
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In the following we prove that the diagonal entries of c+/o<I>(0, (/o+o)) generate analytic semi¬ 
groups when seen as unbounded operators in h 1+ "(§) with dense domain /j 4+a (§), while one of 
the off-diagonal entries is in some sense small (see relation (7.10) below). The analysis is not so 
involved as in the previous sections: on the one hand, the highest derivatives of / and h in (7.8) 
and (7.9) are hidden in the Frechet derivative of the curvature operators (the evolution equation 
(7.4) has a quasilinear structure now), and, on the other hand, we can rely on previous arguments 
and computations. Therefore, we do not prove all our statements in detail. 


The off-diagonal entry. We claim that for each e € (0,1), there exists Kq = A'o(e) > 0 such that 


||+$i( 0 , {f 0 ,h 0 ))\h]\\ 1+a < e 
Recalling (3.12) and observing that 

h" 


4 + 0 ! 


+ A C 


d h K(h 0 )[h\ = 


3 h' 0 ti'ti 


(1 + h’ 2 ) 3/2 (1 + + ) 5/2 

it is useful to split the solution (W+[h], W— [h]) of (7.9) as 

(W+[h], W-[h}) = ( W \, Wi) + c wl , Wl 
where ( IT I, W}_) is the solution to 


i+a for all h € ++“(§). 
for all h £ h 4+ “(§), 


(7.10) 



/ 


+/o> ha)W\ = -d h A(f 0 ,h 0 )[h]v + 

in D + , 





A(f 0 )Wl = 0 

in D_, 


< 

B{fo,ho)W 

l-B(f 0 )Wi = -d k B(fo,h 0 )[h}v + 
Wl-Wi = 0 

on r 0 , 
on r 0 , 

(7.11) 




Wl = gp+h + 3-'f h h , 0 hf ) (l + +) 5 / 2 /i' on Ti, 



< 


Wi = 0 

on r_i, 


while (11+ IF 2 ) solves 






f 

A(f 0 ,ho)Wl = 0 

in n_)_, 





A(f 0 )wl = 0 

in D_, 



< 

B(fo,ho)Wl-B(f 0 )Wl = 0 

Wl-Wl = 0 

on r 0 , 
on r 0 , 

(7.12) 




W 2 = - 7fc (l + h'l)~ 3 / 2 h" 

on Ti, 




k. 

W 2 = 0 

on r_i. 


In view of (3.12) we get 





\\d h $i(0, {fo,ho))[h]\\ 1+a < C( tr 0 VWl pQ + || tr 0 V 1 + || 1+a ) 





<c{\\wl\g- a + \\wl\\£l 2 J n 


(7.13) 

where we again use 

the notation (l/2)fl_ = § x (—1/2,0). Due to (3.6) 

we have 





II Wl\S~ a <C\\h\\ 3+a . 


(7.14) 


Finally, we infer from [1, Theorem 9.3] and the arguments in the proof of Proposition 3.3 that 
\\WltlT- < C(ll^ll„ n+ + IIH^II?-) < C (||^||^ /2 + ||^||^ a/2 ) < C||fc || 1+n/2 (7.15) 

for all h € h 2+a (§>). The relations (7.13)-(7.15) and the interpolation property (4.9) lead us the 
desired estimate (7.10). 
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The first diagonal entry. We start by identifying the “leading order part” dt <h7 € £(/i 4+ “(§), /i 1+ “(§)) 
of 0/$ l( O, (/o, h 0 )). This is defined as ' 

df*l[f]:=-B{f 0 )wl[f], f£h 4+a ( §), 

where {w\ [/], wf_ [/]) denotes the solution to the problem 


•Ao(fo,ho)w+[f] 

= 0 

in 


= 0 

in 

1 

SQ 

1 

^ + 

0 

ST 

= 0 

on Fq. 

wl[f]-wl[f] 

= 7/(1 + 

on Fq. 

<[/] 

= 0 

on Ti, 

wllf] 

= 0 

on r_ 


Here, *4 q(/ 0 ) and A^ifo, ho) are the operators defined by (5.3) (with (/*,/i*) replaced by (fo,ho)). 
The same arguments as in the proof of Lemma 5.1 show that, given e € (0,1), there exists K\ = 
A'i(e) > 0 such that 

\\df^(0, (/„, ho))[f] - df$l[f] || 1+a < e||/||4fa + i^rll/llr+a for all / G /r 4+Q (§). (7.16) 


Theorem 7.2. Lef /i > 0 and a' € (0, a) &e given. Then there exist an integer p >3 , a p-partition 
of unity {nj} 1< / <2 p+i> and a constant = K 2 (p), and for each 1 < j < 2 P+1 there are bounded 
operators A / G £^/i 4+q (§), /i 1+a (§)) swc/i that 


for all f € /i 4+0 


lin^^n/] - A j [n^/]|| 1+a < /xiinf/m+a + ^ll/IU+a' 

The operators A j are defined by the formula 
k 


A j[f] ■= — 


1 + 


/r- V / 0 ~ d 


, Ho [/] - /o(xp tr 0 


where (w+ [/], w_ [/]) denotes the solution to the problem 


.,n,3 


A^(f 0 ,h 0 )w^ [/] 

= 0 

in fl+ 


= 0 

in 

\f\-Bi(fo)wy\f] 

= 0 

on r 0 . 

w^if] - w y\f] 

= v f f" 

on T 0 , 

w^lf] 

= 0 

onT 

w^if] 

= 0 

on r_ 


The operators Aq j(fo,ho), Aoj(fo), 13j(fo,ho), and Bj(fo) are the same as in Theorem 5.5 (with 
(f*,h*) replaced by (f 0 ,h 0 )) and V f := 7/(1 + /o 2 ( 2 y ))~ 3/2 . 

Proof. The proof is similar to that of Theorem 5.5 and is therefore omitted. □ 


From the computations in Section 5 one infers that, given a p-partition of unity {n p } 1< / <2 p+i 
and 1 < j < 2 P+1 , the operator A j is a Fourier multiplier of the form 


A j [ ^2 f me 


AmJmC 

mez 
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the (real) symbol (A m ) m being given by 

A m — — |,(tanh(g + m) + tanh(D.m) N 
* V 02 Dj r m f3 2 D-m ' 

cf. (5.23) and (5.24). Here, and D± are the positive constants defined in Section 5 (with (/*, /i*) 
replaced by (f 0 ,h 0 )). 

Lemma 7.3. Given ( fo,ho ) G Vfl (h 4+ “(§)) 2 and bo G h 2 +a (S), there exist constants > 1 and 
cji > 0 depending only on (fo,ho,bo) such that for any p-partition of unity {n 4 ’} 1<:/< 2 P+i withp > 3, 
the operators A j, 1 < j < 2 P+1 ; satisfy 

A - Aj G Isom(h 4 + “(S), /i 1+ “(§)), 

K l||(A — ^')[/]||l+Q > |A| • ||/||l+a + H/IU+a 

/or all f G /i 4+ "(§) and A G C with Re A > cji . 

Proof. This follows analogously as in the proof of Lemma 5.6. □ 

We are now in the position to prove the generator property for (/o, ho)). 

Theorem 7.4. If (/q, ho) G V n (h 4+a (§)) 2 and bo € h 2 +a (E>), then 

-d f *i( 0 , (/ 0 , h 0 )) G 'H(/i 4+ “(§), /i 1+ “(§)). 

Proof. Because of (7.16) and [2, Theorem 1.1.3.1 (ii)], we only need to prove that — dj&f is an 
analytic generator. The arguments in the proof of Theorem 5.7 together with Remark 5.4, Theorem 
7.2, and Lemma 7.3 show that there exist constants k/ > 0 and uj\ > 1 such that 

«!||(A - d f ^)[f]\\ 1+a >|A| • \\f\\i +a + ||/|| 4 +a (7.17) 

for all A G C with Re A > ui\ and / G h A+a {§). As (A — <9/<3?/) is one-to-one for Re A > ui\ by (7.17), 
we are left to show that (A — dfQJ) is a Fredholm operator of index zero for A > uii, see [2, Remark 
I. 1 . 2 . 1 ], 

We can decompose df&J as df&J = vER o T 2 , where T 2 : h A+a (S) —> h 2 +a (S) is defined by 

T 2 f ■= - — - f" 

(1 + ^ 2 ) 3 / 2/ ’ 

and Tj : h 2 +a (E>) —>• h 1+ “(§) is the operator obtained from df&f 0 [f] when choosing A p = 1 , cf. 
(5.7) and (5.8). A simple consequence of Proposition 5.8 is that is a Fredholm operator of index 
zero, the same Fredholm property being valid also for 'k 2 . By a classical result [34, Theorem 13.1] 
we then know that also df&f = Ti o \k 2 is also a Fredholm operator of index zero. This property 
allows us to conclude that (A — <9/ < 3?/) is bijective for all A > and to finish the proof. □ 

Remark 7.5. Let (/*, h*) G V be given and a G h 2+a ($) be a negative function. Moreover, let A be 
the operator defined by (5.7). Then, 

—A o <9 2 G H(h A+a {S),h 1 +a (S)). 


The second diagonal entry. We define the operator 8^2 ^ £(/i 4+ "(§), /i 1+Q (§)) by the formula 

8 h $m := -B 1 (fo,ho)Wl[h], h G h 4 +“(S), 
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where WT[h] denotes the solution to the problem 

( AUfo,h 0 )W^[h] = 0 mSl+, 

{ Wl\h] =0 on r 0 , 

{ Wl\h] = — 7ft (l + h$)- 3 / 2 h" on h, 

and with *4. 4 (/o,ho) defined by (6.3) (with (/*,h*) replaced by (fo,ho)). The operator dh <&2 is the 
“leading order” part of the Frechet derivative ^^(O, (/o, ho)) in the sense that, given £ € (0,1), 
there exists K 3 = K 3 {e) > 0 such that 

\\dhM 0 , (/o, h 0 ))[h] - d h $Z[h] ||i+a < e||h|| 4+ a + K 3 \\h\\ 1+a for all h € h 4+ “(S). (7.18) 

The estimate (7.18) is obtained by arguing as in Lemma 6.1. 


Theorem 7.6. Let /i > 0 and a' € (0, a) be given. Then there exist an integer p > 3, a p-partition 
of unity {IIj} 1< j < 2 P+i> and a constant iL 4 = K^(p), and for each 1 < j < 2 P+1 there are bounded 
operators Bj € £(h 4+ "(§), h 1+ “(§)) such that 

||n^^[h]-B J [n^]|| l+a L /X||n^|| 4+ a+^ 4 ||h||4+a' 


for all h € h 4+Q (8). The operators B j are defined by the formula 


B j[h\ :=- 


k ( 1 + h'f 


p+\h 0 - /o 


tri d y W+ J [h] - h' 0 (x P j) tri d x W_ 


7rj 



where Wff 3 


[h] denotes the solution to the problem 


' Alj(f 0 ,h 0 )Wl’ j [h] 
< Wl' 3 [h] 

WZ’ 3 [h] 


0 in fl+, 

0 on To, 

-V h h" onTi, 


with AJ j(fo,ho) as in Theorem 6.3 (with (/*,h*) replaced by (fo,ho)) andVh := 7 / l (l + /ig 2 (xp) 3 / 2 . 
Proof. The proof is similar to that of Theorem 6.3 and is therefore omitted. □ 


The computations in Section 6 show that, given a p-partition of unity {nl’} 1<J< 2 P+i and 1 < j < 
2 P+1 , the operator B j is a Fourier multiplier of the form 




= £ (p m h m e imx , 

m$zL 


the (real) symbol (<p m ) m being given by 

„ kVh m 3 

^ tanh(Hm) ’ 

cf. (6.18) and (6.19), where D is the positive constant defined in Section 6 (with (/*,/&*) replaced 
by (fo, ho)). 

Lemma 7.7. Given (/o,ho) € Vfl (h 4+ “(§)) 2 and feo £ h 2+a (§), there exist constants K 2 > 1 and 
u >2 > 0 depending only on (fo,hoAo) such that for any p-partition of unity {rF} 1 < J<2 p+i withp > 3, 
the operators B j, 1 < j < 2 P+1 , satisfy 

A — Bj € Isom(h 4+ “(S), h 1+ “(§)), 


« 2 ||(A - Bj)[/l]||l+a > |A| ■ \\h\\ 1+a + ||h|| 4+a 
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for all h £ h 4+Q (§) and A £ C with Re A > 

Proof. This follows analogously as in the proof of Lemma 6.4. □ 


We now establish the generator property for < 9 /i 4 > 2 ( 0 , (/o, ho)). 

Theorem 7.8. If (fo, ho) £ V H (h 4+Q (§)) 2 and bo € h 2+a (E>), then 

~d h M 0 , (/o, h 0 )) € H(h 4+Q (§), h 1+ “(§)). 

Proof. The proof follows by combining the arguments used to establish Theorems 6.5 and 7.4. □ 


We conclude this section with the proof of Theorem 7.1. 


Proof of Theorem 7.1. In view of Theorems 7.4, 7.8, relation (7.10) and [2, Theorem 1.1.6.1 and 
Remark 1.1.6.2]), we get that 

~d {f , h) m (/o,h 0 ))£^((h 4+ “(S)) 2 ,(h 1+ “(§)) 2 ) 

for all (/o, ho) £ Vn(h 4+ "(§)) 2 and bo = b( 0) £ h 2+ “(§). The remaining part of the proof is identical 
to that of Theorem 2.1. □ 


Appendix A. 

We collect here some explicit formulae for operators used in the previous sections. Given (/, h) € 
V, it follows readily from the definitions in Section 3 of the operators A(f), A(f,h), 13(f), and 
B(f,h ) that 


A(f) =d xx - 2 (1 „ +y)// a 


(l+y) 2 / /2 + l g (1 (f-d)f"-2f» 

\2 u w \ L ^y) ( f j \2 u y> 


f-d ' (f-d) 2 yy v (f-d) 2 

2 1 yh" + (1 - y)f" 


A(f, h) =d xx - 2 Vh ' + } 1 , V)f ' d xy + ^ + . ( * +l d yy - 
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and therefore their (partial) derivatives at a fixed (/*,/i*) € V are given by 
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